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Abstract 

We study a variational model from micromagnetics involving a nonlocal Ginzburg-Landau 
type energy for S 1 -valued vector fields. These vector fields form domain walls, called Neel walls, 
that correspond to one-dimensional transitions between two directions within the unit circle S 1 . 
Due to the nonlocality of the energy, a Neel wall is a two length scale object, comprising a core 
and two logarithmically decaying tails. Our aim is to determine the energy differences leading to 
repulsion or attraction between Neel walls. In contrast to the usual Ginzburg-Landau vortices, 
we obtain a renormalised energy for Neel walls that shows both a tail-tail interaction and a 
core-tail interaction. This is a novel feature for Ginzburg-Landau type energies that entails 
attraction between Neel walls of the same sign and repulsion between Neel walls of opposite 
signs. 

Keywords: Neel walls, Ginzburg-Landau, nonlocal, renormalised energy, interaction, micro¬ 
magnetics 


1 Introduction 

In this article, we analyse a variational model describing the formation of domain walls in ferromag¬ 
netic thin films. These domain walls are called Neel walls and represent one-dimensional transition 
layers connecting two directions of the magnetisation within the unit circle S 1 . Due to dipolar effects, 
the variational problem is strongly nonlocal and generates Neel walls with an interesting core-and- 
tail structure. Our aim is to study the repulsive or attractive interaction between the domain walls 
in terms of their energy. This interaction energy governs the location of the domain walls and is 
analogous to the renormalised energy in Ginzburg-Landau type problems (see the seminal book [3]). 
Although our analysis builds to some extent on the theory of Ginzburg-Landau vortices, our model 
has novel features that have not been studied before. In contrast to the usual Ginzburg-Landau 
vortices, we obtain a renormalised energy for Neel walls incorporating two types of interaction: a 
tail-tail interaction and a core-tail interaction. This is due to the nonlocal character of the model 
and the two distinct length scales of the core and the tails (of logarithmic decay). Moreover, Neel 
walls of opposite signs repel each other and Neel walls of the same sign attract each other, whereas 
Ginzburg-Landau vortices show the opposite behaviour. This observation is consistent with the 
physical prediction (see m Section 3.6.(C)]). Furthermore, in typical Ginzburg-Landau systems, 
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most of the energy is contained in the highest order term, whereas in our model, it is the lowest 
order term that contains most of the energy. From a technical point of view, the lack of a quantised 
Jacobian is an additional difficulty in the analysis of our model. 

1.1 The model 

The magnetisation We consider a one-dimensional model for transition layers (incorporating 
several Neel walls) in the magnetisation of a thin ferromagnetic film. The magnetisation is repre¬ 
sented by a continuous map 

m : (—1,1) —> S 1 . 

More precisely, we can think of a ferromagnetic thin film of the shape (—1,1) x (0, h) x K. (with very 
small thickness h > 0 in the a^-direction) and a magnetisation vector field M : (—1,1) x (0, /i) x M —»• 
§ 2 of the form M( xi,£ 2 , 2 : 3 ) = (m(x i),0). Here, the non-dependence of M on X 2 is a natural 
assumption for a thin film, whereas the non-dependence on x% represents a simplification of the 
problem. (It implies that the walls appear in planes parallel to the X 2 X 3 ~p\ane and we assume that 
the magnetisation depends only on the normal direction aq.) The strip over X\ £ (—1,1) does not 
necessarily represent the whole ferromagnetic sample, but merely a region that contains the Neel 
walls in question. The assumption that the third component M 3 vanishes is consistent with the 
fact that Neel walls correspond to an in-plane magnetisation. Another characteristic feature of Neel 
walls is that the magnetisations on either side (represented by m(— 1 ) and m(l) in our model) differ 
by a vector parallel to the wall plane (in this case the ^-direction). Thus there exists a number 
a £ (0, 7 r) such that 

TOi(— 1) = mi(l) = cosa. (1) 

Moreover, we will sometimes assume that 

m(— 1 ) = m(l) = (cos a, sin a), ( 2 ) 

so that a winding number can be defined. 


m 
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Figure 1: A magnetisation m = (mi,m 2 ) of winding number —1 consisting of a positive Neel wall 
of angle 2a and a negative Neel wall of angle 2 ( 7 r — a) (right). 


More precisely, since m is continuous, there exists a continuous function ip : (—1,1) —»• M, called 
a lifting of m, such that 

m = (cos p, sin ip) in (—1,1) 

and ip(— 1) = a. If © holds, then the winding number (or topological degree) of m is defined as 


deg(m) 


yq)-y(-i) eZi 


The angle a £ (0,7r) will stay fixed throughout this paper. (The case a £ {0,7r} is geometrically 
different and is not studied here.) However, our arguments do not require that TO 2 (—1) = 7712 ( 1 ) = 
sin a in principle and we will present our results in a wider generality, i.e., with 77i2(±l) £ {± sinf?}. 
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The energy The energy for our model comprises two terms, called the exchange energy and the 
magnetostatic energy (or stray-field energy), respectively. The exchange energy is modelled by the 
following expression involving the L 2 -norm of the derivative m'\ 


e 

2 



\m'\ 2 dx i 



e 

2 



Ki ) 2 

1 — m\ 


dx i. 


Here e > 0 is a ratio between a material constant called the exchange length and the length scale 
of the thin film. (This is a model obtained after rescaling, i.e., the length scale of the ferromagnetic 
sample has been set to unit size.) The number e is assumed to be small, and we will eventually 
study the limit e \ 0. 

We write x = (xi, x 2 ) for a generic point in the upper half-plane R+ = lx (0, 00 ). In order to 
compute the magnetostatic energy, we need to solve the boundary value problenij 


A u = 0 
du 
dx 2 
du 
dx 2 


= 0 


in K^_, 

on (-1,1) x {0}, 

on (—oo,—l) x {0} and on (l,oo) x {0}. 


Equivalently, if we extend m ± by the constant cos a onl\(-l,l), then 


Vu-V£dx = / m' 1 ((-,0)dxi for every (eC 0 °°( 


(3) 

(4) 

(5) 

( 6 ) 


Let W 1,2 (K 2 ) be the completion of Co°(K 2 ) with respect to the norm 


IICIIvb 1 ’ 2 (R 2 ) — I|VCIIl 2 (M 2 ). 

(We sometimes abuse notation and treat elements of IT 1,2 (M 2 ) as functions, even though the com¬ 
pletion process identifies any two functions that differ by a constant.) For an open set C R 2 , we 
write W 1,2 (fi) for the set of all restrictions of functions in W 1 ’ 2 (]R 2 ) to H and 


IICIIivi.2(n) — II^CIU 2 (n)- 


By the Lax-Milgram theorem, solutions of © are unique in W 1,2 (K^) (i.e., up to a constant). Thus 
the quantity 


1 

2 


|Vw| 2 dx 


depends only on mi. This is the term representing the magnetostatic energy. It is worth remarking 
that the solutions u of © © in TT 1,2 (R^_) have a limit for |x| —» 00 . Indeed, if we extend u to R 2 
by even reflection, then we obtain a harmonic function near 00 with finite Dirichlet energy, and it is 
well-known that the limit exists at 00 . Then we normalise this constant and define U(m) (sometimes 
also denoted U(m 1 )) to be the unique solution of © in IT 1,2 (R^_) with 


U(m) —> 0 as |x| —> 00 . 


Moreover, in view of ©, using the extension of m\ by the constant cos a onR\(—1,1), we may 
express the magnetostatic energy in terms of the homogeneous || • ||^i/ 2 -seminorm of mi (see e.g. 

BUM)- 


1 

2 


f „ 1 r 

d 

1/2 

/ |VC/(m)| 2 dx = — 

/ r 2 . z Jr 

dx 1 

m\ 


dx\. 


(7) 


1 Here. V u represents the stray-field associated to M, which is also invariant in the 23-direction. 
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To summarise, we study the energy functional 

E e {m) = ^ f \m'\ 2 dx\ + 7 - f \S7U{m)\ 2 dx 

2 J -1 2 Jri 

for to £ W 1,2 ((—1, l),^ 1 ) satisfying (JT]). We are interested in the behaviour of to and of its energy 
E e {in) as e \ 0, especially under conditions that force the nucleation of several Neel walls. 

Neel walls If we trace to from —1 to 1, we may well find that m winds around the circle S 1 one 
or several times. If that happens, then there necessarily exist two points a+,a_ £ (—1,1) such that 
77Ti(a+) = 1 and mi(a_) = —1. But even if the topology of to is trivial (i.e., if deg(m) = 0), a 
transition from (cos a, sin a) to (cos a, — sin a) may occur, giving rise to a point in between where 
7711 reaches one of the values ±1. We think of any such transition as a Neel wall and we use these 
points in order to track them. Obviously, it is possible for toi to attain ±1 when no proper transition 
occurs, but from the energetics point of view, this makes no difference and we call this a Neel wall 
anyway. We speak of a positive or negative Neel wall depending on the sign of toi (see Figure [[]). 
We will see that a Neel wall has a two-length scale structure comprising a core of size 5 = elog - 
around the transition point and two tails of size 0(1), where toi decays logarithmically to cos a 
(see Theorem [55] below). The total change of the phase during the transition is called the rotation 
angle of the Neel wall (which may be 0 by the above convention) @ For more physical background, 
we refer to nnnni- 

We will assume in the following that there are certain points ai,..., ajv £ (—1,1) such that 


— 1< «!<•••< ajv<l (8) 

and certain numbers d\,, djy £ {—1,1} such that 

mi (a n ) = d n for n = l,...,N. (9) 

These points (a ra )i<n<jv represent the positions of the Neel walls that we study, while (d ra )i<n<iv 
indicate whether a Neel wall is positive or negative. We keep the number N of walls fixed throughout 
the paper. Let 

An = {a = (asi,..., ajv) £ (—1,1)^ with ©} . 

For a £ An and d £ {±1}^, we consider the set 

M(a, d) = {to £ W 1 ,2 ((—1,1); S 1 ) with (CQ) and {9])} . 

Our aim is to answer the following question. 

Question. For a given a £ An and d £ {±1}^, what is the behaviour of 

inf E e as e \ 0? 

M(a,d ) 

That is, if we prescribe Neel walls at the positions 01 ,..., on with signs d\,. .., dN, what energy 
does it take to achieve such a configuration? We first note that a minimal configuration to always 
exists and that its first component toi is unique. (Obviously, |to 2 | is also unique, but the sign of 
the TO 2 component can change between a„ and a „+1 for two different minimisers.) 

2 When studying the interaction between a pair of walls, the physics literature (see m) distinguishes between 
winding walls, which refers to a pair with the same rotation sense, and unwinding walls, which refers to a pair 
with opposite rotation sense. Except for degenerate cases, a pair of Neel walls with opposite signs according to our 
terminology corresponds to winding walls and a pair with the same sign corresponds to unwinding walls. 
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Proposition 1. There exists a minimiser of inf m (a, d.) E e for any e > 0. Moreover, any minimiser 
m is smooth on (—1,1) \ {ai,..., on} and, has a unique m\-component. 

Proof. The direct method in the calculus of variations yields a minimiser m of E e in M(a,d). The 
regularity of m is standard (see, e.g., [14]). The uniqueness of nil follows from the strict convexity 
of 0 and of the function ( v , w ) >->• 1 f_ w2 for ( v , w) S ffi. x (— 1,1). □ 

We look for an expansion of infj\/(a,d) similar to [9], where it is shown that 


N 


inf 

M(a,d) 


^ = £ 

n= 1 


7 r {d n — cos a) 2 


2 log 7 


+ 0 


( log log 7 \ 

IN \f) 


( 10 ) 


for e > 0 small. Since this is not good enough to understand the interaction between domain walls, 
we need to determine the second term in such an expansion completely and identify the third term 
as well. This problem is analogous to finding the “renormalised energy” in Ginzburg-Landau type 
problems, but in the context of Neel walls, it has remained open until now. 

We give the answer to this question in Theorem [2j The key is to identify the contributions to the 
renormalised energy coming from the interaction between two tails and between a core and a tail of 
two different walls. It turns out that the above expansion is easier to understand when we replace 
e by 6 = elogi (recall that this is the typical length scale of the core of a Neel wall). The first 
two terms of the expansion (U0[) are then united in a single leading order term in the expansion in 
1/| log<5| (see dill) below). The next-to-leading order term corresponds to the renormalised energy. 


1.2 Motivation 

There are several reasons for asking the above question. First, we may want to study the positions of 
Neel walls in equilibrium. Once we have determined the renormalised energy, we can find the likely 
positions by minimising it. Second, we may want to study the dynamics of Neel walls (see, e.g., 
mi). The dynamics of the magnetisation is described by the Landau-Lifshitz-Gilbert equation, 
which is derived from the micromagnetic energy through a variational principle. For reasons that 
are explained below, understanding the asymptotic behaviour of the energy is expected to be a 
key step towards deriving an effective motion law for the walls in the limit e 0 . A third reason 
for studying these long range interactions is that we want to understand some phenomena in thin 
ferromagnetic films where they matter, such as cross-tie walls. A cross-tie wall is a typical domain 
wall that consists in an ensemble of Neel walls and micromagnetic vortices (similar to Ginzburg- 
Landau vortices), see 0 CD J2E iiSJ- It has an internal length scale, the size of which is not predicted 
by any existing theory, and our analysis on the interaction energy of Neel walls could represent an 
significant step forward here. 

A related question concerns the analysis of general transition layers m carrying a winding number 
when the location of the Neel walls is no longer prescribed. More precisely, suppose that the 
lifting q) : (— 1 , 1 ) —> R of m = (cos </?, sin ip) satisfies the boundary conditions ip(— 1 ) = a and 
</?(l) = 2 £- 7 t + a, so that we have winding number i , i.e. t = deg(?n). Hence the magnetisation 
performs i full rotations, so that 0 is satisfied. Then by continuity, we necessarily have a certain 
number of transitions between (cos a, sin a) and (cos a, — sin a). 

Open problem. For a prescribed winding number and given suitable control of E e (m), what can 
we say about the profile of m and of the stray field potential U ( m ) ? 

As mentioned before, a prescribed degree £ will automatically give rise to certain Neel walls. But 
it is not obvious, for example, that these Neel walls stay separate from one another (uniformly as 
e —> 0) and that one can rule out other transitions. In fact, it is an open question whether the lifting 
of to is monotone even for minimisers (which would exclude unexpected transitions). However, 
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assuming good control of the energy, we expect to have exactly 2 £ transitions (corresponding to the 
expected Neel walls) and no extraordinary behaviour of the magnetisation in between. For the stray 
field energy, it is expected that the energy density concentrates at the walls. Such information would 
be useful in the study of compactness properties in the appropriate function spaces, for example 
with a view to T-convergence. 


1.3 Main results 


For any e £ (0, |], let 

S = e log - 
e 

and define the metric g on (— 1 , 1 ) bjH 

e(b, c ) = g [0) 1) for b,c £ (—1,1). 

We have the following result, answering the question on page [I] 

Theorem 2. There exists a function e : {±1} — > R such that for any a £ An and d £ {±1}^, the 
following holds true. Let 7 n = d n — cos a for n = 1,..., N and let 


N 

£■ 

n— 1 


(a,d) = } e(d n ) - 


7r 


N 

71=1 


N 


log (2 - 2 a 2 n ) - - 7kln log 


n= 1 k=fin 


1 ~b \/lp(dfc, Un)^ 
£»(afc, o n ) 


Then 


W(a,d) = lim 

e \,0 



„ 7T 
inf E e — — 
M(a,d) 2 



In analogy to the theory of Ginzburg-Landau vortices, we call W (a,d) the renormalised energy 
for the N walls placed at a = (ai,...ajv) with signs d = (di ,..., d/v). As the theorem shows, 
W(a, d) represents the next-to-leading order term in the expansion of infM(a,d) E e in 1/| log<5|. If we 
express these asymptotics in terms of e, our result improves m by determining the precise second 
and third coefficients: 


inf = ~ 1 2 (l-J^Tridn-cosa) 2 ("log - + loglog +W(a,d)) +o( 1 ) . (11) 

M(a ’ d) (fog 7 ) \ 2 n=i V e e 2 J \ (fog 7 ) J 

We now briefly discuss how the above expression comes about. Suppose that for a given a £ An, 
we study minimisers m of E t in M(a,d). When e is small, we expect to have a typical Neel wall 
profile near each of the points ai,..., ojv with the prescribed signs di,, dN, and the full transition 
layer m is essentially a superposition of all of these. As discussed previously, we can think of a Neel 
wall as consisting of two parts: a small core around a„ and two logarithmically decaying tails. In our 
situation, the walls are confined in the relatively short interval (— 1 , 1 ) and each tail will interact with 
the other walls and with the boundary as well. We can then account for the full energy inf M(a,d) E e 
(at leading and next-to-leading order) as follows. 


Core energy. The core of each wall requires a certain amount of energy, namely 

e(l) J e(-l) 

-o and - 7T 

(log i) (log i) 

3 We will not use the fact that q is a metric, but if we want to verify it, we can use that £(3^(6), ^(c)) = g(b, c) 
for the Mobius transforms defined in (1261) below for every d E (—1,1)- For the triangle inequality, it then suffices 
to show that g(c , 0) < g(b , 0) + g(b , c) for b, c E (—1,1), which is not difficult. 
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for a positive and a negative wall, respectively. The constants e(±l) are given in Definition 1261 
below as limits of a rescaled energy of the core profile as e —> 0. Then the sum accounts for 
the term 

E»=l <dn) ' 

( lo Sl ) 2 

This is the only term where we have a contribution from the exchange energy and it appears 
only at next-to-leading order in the full energy. All the remaining terms below come from the 
magnetostatic energy alone. 

Tail energy. The two tails of the wall at a n give rise to the energy 

2 log*’ 


leading to a total of 

2 log* 

This is the leading order term of the full energy. 

Tail-boundary interaction. Moving a wall relative to the boundary points ±1 will deform the 
tail profile, resulting in a change of the energy. This phenomenon gives rise to the energy 


Tfo l l°g(2 ~ 2 a 2 n ) 

2 N |) 2 


for the wall at a n . Summing up these contributions, we obtain 


7r 

2 (log l ) 2 


N 

Xb 2 l° g (2- 2<4)- 

n=1 


(The sign here is not a mistake; it is the opposite of the sign of the corresponding expression 
in Theorem [2]) This means that the tails are attracted by the boundary, in the sense that the 
energy decreases if a n approaches ±1. 

Tail-tail interaction. There is an energy contribution coming from reinforcement or cancellation 
between the stray fields generated by different walls. For the walls at ak and a n with k ^ n, 
this amounts to 

TTJkJn j | 1 + \A ~ g(«fc,Qrt) 2 

2 (fog*) 2 \ Biflkidn) 

The total contribution is 



7r 


N 


2 (log*) n =l fc/n 


log 


1 + -y/1 - e(afc,«ra) 2 
Qi'Q’ki ^n) 


(Again we have the opposite sign relative to the above theorem.) A conclusion is that the tails 
of two walls attract each other if they have opposite signs and repel each other if they have 
the same sign0 

4 This is because the function t i-» - is decreasing on (0,1). 
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Tail-core interaction. Since the profile of a Neel wall decays only logarithmically, it will change 
the turning angle of the neighbouring walls slightly. This has an effect on the energy as well 
(at the next-to-leading order). Indeed, the tail of the wall at a,k and the core of the wall at a n 
with k ^ n lead to a contribution of 

TT7fc7n lo ( 1 + Vl - g (q fc ,q,) 2 \ 

(log |) 2 \ g(ak,dn) J 

We also have an interaction between the two tails of a wall and its own core: if k = n, then 
we obtain the energy 

*1 2 log(2 ~ 2a 2 ) 

(log l ) 2 

This gives a total of 


N 


(log l) n=l 


7 2 l°g(2 - 2a 2 ) 


k^n 


Ikln log 


1 + y/1 — g(Qfc, Qra) 2 

g{ak,a n ) 


This is twice the size of the terms from the tail-boundary interaction and tail-tail interaction, 
but with the opposite signs, resulting in a net repulsion between walls of opposite signs and 
a net attraction between walls of the same sign. Furthermore, we have a net repulsion of the 
walls by the boundary. 


Notwithstanding the term ‘energy’ used in this description, strictly speaking, these are energy 
differences and therefore some of them may be negative. All except one of these contributions occur 
similarly in the theory of Ginzburg-Landau vortices. The core-tail interaction, on the other hand, 
is new and more delicate to handle. 


I. 4 Physical relevance 

Our result represents a rigorous proof of the physical prediction on the interaction energy between 
Neel walls. Indeed, Hubert and Schafer m Section 3.6. (C)]) predict the following behaviour 
in the case of a pair of Neel walls: “The extended tails of Neel walls lead to strong interactions 
between them [...] The interactions become important as soon as the tail regions overlap. The sign 
of the interaction depends on the wall rotation sense. Neel walls of opposite rotation sense (so-called 
unwinding walls) attract each other because they generate opposite charges in their overlapping tails. 
If they are not pinned, they can annihilate. Neel walls of equal rotation sense (winding walls) repel 
each other. ” (We recall that unwinding walls correspond—according to our definition in Section 

II . 11 -to a pair of Neel walls with the same sign, while winding walls correspond to a pair of walls 
with the opposite signs as in Figure ID) 


1.5 Comparison with a linear model 

If we replace the exchange energy by the simpler expression 


e 

2 



dx i, 


then the energy functional, now given by 


E e (mi) 


(m[) 2 dx i 




|VI/(?ni)| 2 dx, 


nii '■ (— 1,1) — > K, 









becomes a quadratic form and the Euler-Lagrange equation for its critical points becomes linear. 
This functional has been used as a tool for studying the energy of Neel walls mm- Since the 
exchange energy in E e does not enter the expansion m at the leading order, we may expect good 
approximations from the linear model involving the energy functional E e . The exchange energy does 
have an effect on the next-to-leading order term however, even though it is not through a direct 
contribution but rather by changing the core width of a domain wall. For the linear model, the 
core width of a domain wall is of order e. Accordingly, for the functional E e , the expansion that 
corresponds to (flT)ll is of the form 


inf E e (mi) = 

m£M(a,d) 


N 

E 

n= 1 


Tr(d n — cos a ) 2 W(a,d) 


2 log 7 


(l°g 7)^ 


. ( lo S e ) ‘ 


as e \ 0. Here, W is nearly the same as the function W from Theorem [2j except that it may differ 
by a number depending only on N and d. That is, there exists a function e : {±1} —> R such that 


(a,d) = 


N 

E 

n= 1 


e(d n ) 


N 

- -E 

2 ^ 

n= 1 


N 


7 1 l°g (2 - 2 a 2 ) - -^E 7fc7n log 


n= 1 k^n 


1 + y/l- gjak^n) 2 

diP'k , a n ) 


As for the full model, we may regard e(±l) as the core energy of a transition of sign ±1. Our 
analysis does not give an explicit expression, but for variational principles where the number and 
signs of the Neel walls does not change, this part of the limiting energy is irrelevant. 

The formula can be proved with the same arguments as in the proof of Theorem[2]below, although 
the linearity allows a few shortcuts. Therefore, we do not give a separate proof but leave it to the 
reader to make the necessary changes. 

As a consequence, the linear model does not describe the interaction between Neel walls accu¬ 
rately, but the discrepancy is easily corrected by adjusting the core width (i.e., replacing e by <5). 
Although we study only the energy of interacting Neel walls in this paper, the analogy to the theory 
of Ginzburg-Landau vortices (see Sect. 11.61) suggests that the same may be true for the dynamics 
of Neel walls. The simplified model may therefore be useful as a test case for future analysis, or, 
with the necessary care and the appropriate corrections, even be used for quantitative predictions. 


1.6 Comparison with Ginzburg-Landau vortices 


The interaction between topological singularities has been intensively studied in the last two decades 
in the context of Ginzburg-Landau problems. The work was pioneered by Bethuel, Brezis, and 
Helein mm, and an overview of later developments can be found in a book by Sandier and Serfaty 
[31] . These problems are designed to describe phenomena in superconductors and Bose-Einstein 
condensates, and a simple model that captures some of the main features is based on the functionals 

Gt(/) = J n Q |v/|2 + h {l ~ l/|2)2 ) dx (12) 

for a domain Slcl 2 and a function / : Q —> R 2 . We identify R 2 with C. Then in the limit e \ 0, 
the analysis in the aforementioned papers leads to a limiting function / : Q —> C of the form 


/(*) 



d n 


for certain points a ±,... ,un 6 B, integers di ,... ,dw € ^\{0}, and a function 9 : —> R. The 

renormalised energy 


lim inf 

r\0 



'^\Un=l B rM 


IV / 1 2 dx 
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appears in a result similar to Theorem [2] (together with an additional term describing the core en¬ 
ergy). Here B r (a) stands for the open ball of radius r and centre a. We have topological singularities 
at ai,..., a n with vortex structures and with topological degrees d \,..., d n . These data are encoded 
in the distributional Jacobian 

J(/) = icurl(/ J --V/), 

where f ± = (-f 2 ,fi)- 

The renormalised energy gives information about the vortex positions in equilibrium, but it is 
also important for their dynamics. Typically, if / evolves by a variational equation derived from G e , 
then on an appropriate time scale, the limiting motion law for the vortices (as e \ 0) is described 
by an analogous equation derived from the renormalised energy. This is true for gradient flows 
ED [Ml HU, Schrodinger type equations 0E1, as well as nonlinear wave equations Ei mi- 

it has been observed before that certain phenomena from micromagnetics give rise to similar 
models DU ng ng ini Eg. The connection to our model is less obvious, but can be seen once we 
show that under assumptions such as in Theorem H we obtain a limiting function from the rescaled 
stray field potential (log j)U(m ) of the form 

<dO) = M x ) + ( arctan (———) - 

“ V \xi-a n J 2\xi — a n \ 


for some a € An and d e {±1} W and a harmonic function it* : Mi —> M that is smooth near 
(—1,1) x {0} (see Sect. [2] for details). Examining u* d near the point (a n , 0) € M 2 , we see that it 
behaves like the phase of a vortex in the upper half-plane, up to the constant 7 „. The expression 


hmin f ( \S7u* d \ 2 dx 

r V> \ Z JR 2 + \\J% =1 B r (a n ,0) 



also plays a role, although for our problem, it only accounts for a part of renormalised energy in 
Theorem [2] (even after adding the core energy). In fact, a Neel wall at a n behaves like a vortex 
of “degree” 7 „ in many respects, which is why the toolbox from the theory of Ginzburg-Landau 
vortices is very useful for the analysis. 

There are, however, significant differences to Ginzburg-Landau vortices as well. A Neel wall is a 
two-length scale object, comprising a core and two tails, each with its own characteristic length. In 
contrast, in the standard Ginzburg-Landau problem, a vortex has a single length scale characterising 
its core and the renormalised energy between the vortices comes essentially from the interaction of 
its out-of-core structure. For Neel walls, we have a renormalised energy consisting of two parts. The 
interaction between the tails of two walls is similar to the interaction between Ginzburg-Landau 
vortices and gives rise to the above expression. But in addition, we have an interaction between the 
core of one wall and the tail of another, which is a novel feature for Ginzburg-Landau type systems. 
This interaction is responsible for the fact that we have attraction for walls of the same sign and 
repulsion in the case of opposite signs, whereas for Ginzburg-Landau vortices, we have attraction 
for degrees of opposite signs and repulsion for degrees of the same sign. Finally, our “degree” 7 „ 
is not quantised in the same way as the degree of Ginzburg-Landau vortices. It does take only 
two values (±1 — cos a), but these depend on the choice of the angle a and are not topological 
invariants. As a consequence, the Jacobian becomes a much less powerful tool. To overcome these 
difficulties, we use duality arguments and “logarithmically failing” interpolation inequalities (see 
musi), T-convergence methods, and refined elliptic estimates. 

In our model, the magnetostatic energy, being of order 0(1/| log <5|), dominates the higher order 
exchange energy (of order 0(l/(log<5) 2 ) for small values of e. This is in contrast to most Ginzburg- 
Landau systems, where the highest order term is dominant. This is the case for the functionals 
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G e in (1121) . but also for similar models coming from micromagnetics. For example, the model for 
boundary vortices studied by Kurzke |18] contains a term coming from the exchange energy and one 
coming from the magnetostatic energy as well, but in the analysis, the roles of the two are reversed 
relative to the model for Neel walls. 


1.7 Notation 

We now introduce some notation that we will use frequently throughout this paper. As mentioned 
previously, we define 5 = elog This is the scale of the Neel walls’ typical core width. 

If xo £ R 2 and r > 0, then we write B r {x o) for the open ball in R 2 of radius r centred at 
xq■ Furthermore, we write B+(xq) = B r (x o) HR+. For a set 5 C R 2 , we also use the notation 
d+S = dS n R 2 .. 

For a vector £ = (£i, £ 2 ) £ R 2 , we write £ x = (—£ 2 , £ 1 ). If / : fl —> R is a function on a domain 
n C R 2 with gradient V/ = (|£, J£), then V^/ = (-|£, J£). 

For a £ An, we define 

p(a) = — min{2ai + 2, a 2 — ai,..., ajv — ajv-i, 2 — 2a«}. 

Thus this is a quantity that controls the distance between two points of a and the distance to the 
boundary. For r > 0, also define 

N 

B*M= (J B+(a n , 0) 

n =1 

and 

n r (a) = R 2 \l?*(a). 

Given a function / : R+ —> R with a well-defined trace on (—1,1) x {0}, we often use the 
shorthand notation /' for and 

J f dxi = j /(xi,0) dx 1. 


In addition to the space IT 1 , 2 (R 2 ) introduced in Sect. 11.11 we also define bF*’ 2 (R^; a) for a £ An, 
which is the space of all u £ fUo^OMa)) such that 


l|V«|li a( n r (a)) 

sup--V 11 

r >0 log (- + 2 ) 


< OO. 


If b £ (—1,1) is a single point, then W* 1 , 2 (R^; b) is defined similarly. 


2 The renormalised energy 

When we study minimisers of E e in M(a, d ) and let e tend to 0, then we expect the suitably rescaled 
stray held potential to converge to a harmonic function with specific boundary data depending on 
a and d. We compute this function here, in order to obtain the limiting magnetostatic energy. This 
corresponds to the sum of the tail self-energies and the contributions of the core-tail, tail-tail, and 
tail-boundary interactions. As a side product, we will also obtain information about the expected 
logarithmic profile of the Neel walls. 
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2.1 The limiting rescaled stray-field potential 

Fix a £ Ajsi and d £ {iljW. We recall that 


7 n = d n — cos a for n = 1,..., N 


N 


and we denote 

= f lk ~ ) ’ n = 0,...,N. 

\k=n+l k =1 / 

We look for a solution of the following boundary value problem: 


A< d = 0 


*■+> 


du 


i* a ,d = cr o on (-l,ai) x {0}, 

l *a,d = a n On (a n , a n+1 ) for n = 1,. .., IV - 1, 

L *a,d = a N On (a N , 1) x {0}, 

* 

— = 0 on (—oo, —1) x {0} and on (1, oo) x {0}. 


8X2 


(13) 

(14) 

(15) 

(16) 

(17) 


The boundary data are chosen so that we have a jump of size — 7 T 7 „ at a n for n = 1,..., N. We also 
require that u* d £ Wl ’ 2 (R 2 -a) so that the boundary conditions make sense. However, the Dirichlet 
energy cannot be finite on all of since any solution u* d will behave similarly to the functior 0 


N 


y 7 n I arctan 


X2 


Xi - a r 


7 r(a;i - a n ) 
2\x± a n | 


near the boundary (—1,1) x {0} (see Proposition [5] below). But we can compensate by considering 
the expression 


K,d( u *a,d ) = \ iiminf [ f |V< d \ 2 dx - Trlog - y 7 ^ ) 

2 -NO yj Qr(a) J 


(18) 


and this will be part of the limiting energy. We will see that the unique solution it* d of m mi) 
is the minimizer of E* d and corresponds to the limit as e \ 0 of the rescaled stray field potential 
(log j)U(m e ) associated to the magnetisation m t with walls at ai,...,ajv with prescribed signs 
g?i, ..., d/v (see Proposition [5] and Theorem [22]) . 

In order to solve PD El' we first study the simpler problem 


Ait = 0 
7r 




u = y on (- 1 , 0 ) x { 0 }, 
u = ~\ °n ( 0 , 1 ) x { 0 }, 


du 

—— =0 on (— 00 , — 1 ) x { 0 } and ( 1 , oo) x { 0 }. 

OX 2 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 


We can obtain a solution u* d of El El by a linear combination of solutions to problems of the 
type (fl9l) (l22l) (see Proposition 0. 


5 This function satisfies OH lfl 6 li . but not GB- It does not belong to W t 


i/ 1 ’ 2 (p2 


; a) due to its behaviour at oo. 
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We first construct an explicit solution, using the fact that harmonic functions remain harmonic 
upon precomposition with a conformal map. We identify R 2 with C. Consider the following domain 
in the complex plane C: 


S = {wi + iw 2 £ C: w\ >0 and 0 < w 2 < 7 r} . 
Also consider the conformal map F : S —> R^ with 

1 


F(w) = —- 


w £ S. 


cosh w 

Extend F continuously to the boundary dS \ {^f}. Assuming that u solves (fTill) 


set 


u = u o F. 


Then u solves the boundary value problem 


A u = 0 in S', 

7T 


U = 


„ 7T 
U ~ 


on ( 0 , oo), 
on ni + ( 0 , oo), 


du 

—— =0 on j( 0 , 7 r). 

UX i 


This problem has an obvious solution, 


u(w) = -Im w, 

K ’ 2 


with Imw = w 2 for w = w\ + iw 2 . Thus we obtain a solution of (1191) (1221) by 


u{z) =u(F~ 1 {z)) = - -Imr 1 ^), z e R+, 


which can be written as 


7TZ 

u = Im/ for /(z) = —— F - 1 (z), z £ R^_, 


(23) 


where / is an holomorphic function. Since linrni^oo F 1 (z) = ^r, we conclude that limui^oo u(x ) = 

0 . 

Proposition 3. The function u from (1231) is the unique solution of m ca in W*’2(R^ ;0 ). It 
satisfies |u| < ^ in R+ and limui^oo u{x) = 0. Moreover, it is odd in x\, that is, u(x i,x 2 ) = 
—u{—x 1 ,^ 2 ) for every x € R+. Furthermore, there exists a constant C such that 


Vu(x) 



< C\x\ 


(24) 


for all x £ By 2 {0 ). 

Proof. It is clear from the construction that |it| < ^ and liiri|^ 1 u(x) = 0. Since the Dirichlet 

energy is conformally invariant, we have 



\\7u\ 2 dz 



(Rl\B r (0)) 


|Vd | 2 dw. 
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Note that 


| cosh «;| 2 = - (e w + e~ w ) (eP + e - *) 


2 Re it; _j_ ^—2 Re w _|_ ^ 2ilmw _|_ ^—2ilmw\ 


Thus 


It follows that 


= 4 (e 

= - cosh(2 Re w) + - cos(2 Im w) for w £ C. 
F _ 1 (i?^_\I? r (0)) C 6 S: Re w < ^ arcosh + 1^ | 


(25) 


/ 

JR 2 


|Vw | 2 dz < — arcosh ( — + 1) . 


/M 2 \B r (0) 

In particular, we have u £ Wl ’ 2 (R 2 ; 0 ). 

If u £ W» 1 , 2 (R+;0) is another solution of (fTUl) (121 . then the difference v = u — u belongs to 
B* 1 , 2 (R+;0) and thus to W 1,P (B^ (0)) (~l fT 1 , 2 (M+\-Bi/ 2 ( 0 )) for any p £ [1,2). Furthermore, it is 
harmonic in R^_ and satisfies v = 0 on (—1,1) x { 0 } and = 0 on (—oo,— 1) x { 0 } and on 
(l,oo) x {0}. Standard regularity theory then implies that v £ 4F 1 , 2 (Rj_), and with an integration 
by parts, we obtain 

f |Vu | 2 dx = 0. 

J M 2 

Hence u is the unique solution of CU) (El in W* 1 , 2 (R^; 0 ). 

The odd symmetry of u is a consequence of the uniqueness, as i 4 — u(—x 1 , 2 : 2 ) is another 
solution of the boundary value problem. 

Finally, we consider the function 


v(x) = u( x) — arctan 


X2 


TTX\ 


Xi) 2\xi\ 


which is harmonic in as well with v = 0 on (—1,1) x {0}. Invoking standard regularity theory 
again, we conclude that v £ C°°(By 2 ( 0)). By the odd symmetry, we have Vw(0) = 0. Hence we 
obtain inequality (HO). 

Definition 4. For every b £ (—1,1), we introduce the Mobius transform 


•‘(',-TTE- 

and its inverse We also define Ub : R+ —> R by 

Ub = u o in R^_, 

which, by Proposition 0 is the unique solution of the boundary value problem 


□ 

(26) 

(27) 


A Ub = 0 
7r 

u b = — 


m . 




on (— 1 , b) x { 0 }, 


in the space Wf 




u b = on ( 6 , 1 ) x { 0 }, 
du 

—— =0 on (— 00 , — 1 ) x { 0 } and on ( 1 , 00 ) x { 0 }, 

ox 2 

.;*)■ 
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Now we can also construct a solution of O HU by superposition. 
Proposition 5. The function u* d , defined by 


N 

K,d = J2'ynUa n , 

n =1 


( 28 ) 


is the unique minimiser of E^ d among all u £ W*’ (R+;u) satisfying uni) dm). Moreover, u* d is 
the unique solution of (fl3l) (fI71) inWl’ 2 (M? + \a). 

Proof. By Proposition [U we know that u* ad £ W* 1,2 (R+;a), and it satisfies (fl4l) (fTbl) . Let u be 
another function with these properties. Since it* d is harmonic in integration by parts leads to 



Vu* d 2 dx = f \S7u\ 2 dx — 

[ Vu* d | 2 dx — 2 

[ (Vu — Vu* ad ) 

■ V u* ad dx 

Jn r (a) 

J Q r (a) 

/fi r (a) 


= ( \X7u\ 2 dx — 

[ \Vu* a>d \ 2 dx-2 

[ (u- U* d ) 1 / - 

V< d da. 

J Q, r (a) 

J Q r (a) 

ldn r (a) 



Using (HHl) . we find that there exists a constant C such that for any n = 1 and for x £ 


d K,d , , _ 7n^2 

dxi X (xx - a n ) 2 + x\ 


du *a,d, S Tn(Xl-a n ) 

r -, (*£) / \o o 

OX 2 (Xi-anY + Xi, 


Hence the boundary integral will tend to 0 when we let r 0. Therefore, 


K,d«d) = K,d(u) - tj lim 

Z r\j 0 


Vu- Vu* 


dx. 


Q r (a) 


The limit on the right hand side exists, because the quantity is monotone in r. 

This implies that it* d is the unique minimizer of E* d . The uniqueness of solutions of HU HU 
follows with the same arguments as for Proposition [3] □ 


Remark 6. Since we have an explicit representation of it* d , an easy computation shows that the 
liminf in the definition of E* d (u* d ) is in fact a limit. By the preceding computation, the same holds 
for E* d (u) for any u £ W r * 1,2 (R^_;a) satisfying (fT-H) HU- 


2.2 The energy of the limiting rescaled stray field 

Next we want to compute the energy E* d (u* d ) defined in HU- Since this depends only on a and 
d , we use the abbreviation 

Wx{a,d) = E* a4 {u* a4 ). 

This quantity corresponds to the tail-boundary and tail-tail interaction energy. 

Proposition 7. If g is the metric on (—1,1) defined in Section ] 1.31 then 


Wi{a,d) = 


N 

-Y 
2 ^ 

71=1 


7 n lo §( 2 - 2a l) + 


2 ^ 

k^n 


TkTn log 


1 + y/1 - g(«fc,«n) 2 

g(ak,a n ) 


(29) 


For the proof, we need the following two lemmas. First we compute the rescaled tail-tail inter¬ 
action energy for two Neel walls located at two points b ^ c. 
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Lemma 8. For all b,c£ (—1,1) with b ^ c, if Ub and u c are defined by (1271) . then 

I 

J 


VtJt, • \7u c dx = 7rlog 


1 + s]\ - g{b,c) 2 


g{b,c) 


( 30 ) 


Then we compute the rescaled tail self-energy together with the tail-boundary interaction of a 
Neel wall located at a point b £ (—1,1). 


Lemma 9. There exists a constant C > 0 such that for every b £ (—1,1) and r £ (0,1 — |&|), 


l\Br(b) 


|Vrt(,| 2 dx — 7rlog 


2-2 b 2 


< C 


|6|r r 2 

1 — b 2 (1 — b 2 — \b\r) 2 J 


(31) 


In particular, 


lim 

r\,0 


( f |Vzih| 2 dx — 7rlog - ) = 7rlog(2 — 2b 2 ). 

\jRl\B r (b) rj 


Proof of Lemma\Q Using the definition of the Mobius transform it is easy to check that 

Set q = jE§^- Then |g| = g(b,c) and 


o $ c = <!> b+c . 


/ Vuh-Vu c dx= / Vzt • V(zt c o ${,) dx 
Jr 2 , Jr 2 


Vm • Vu q dx 


i: 


= -PV / U q - — dx\ 


du 

'd ^ 2 


r q du 


nl du 


= -7 T p V/ T^—dx 1 + -PV/ —— dx\. 


l-i dx ; 


dx 2 


In order to determine we recall that the function / defined in (E3l) is holomorphic in R+. Hence 
we have 

... . du du n 

/(2) = te +! te' 




In particular, 


^-=R ef'(z), 2<ER+. 

OX 2 


Differentiating both sides of the equation 


— ~w = f(F(w)), 


we calculate 


f'(F(w)) = —- 


cosh w 


F'(w) sinhw 

Let t £ (0,1) and s = arcosh j. Set w = s + iir, so that F(w) = t. Then coshw = —I and 
sinh w = — sinh s = — \Jt~ 2 — 1. Hence 

/,<i) = wh=? = s ( log |f| “ log 0 + ' /rr75 )) - 
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That is, 


J^(zi,0) = ( log l ;Cl l - log (i + v^ 1 - - *i€(0,l). (32) 

By the odd symmetry of u in x\ (see Proposition [3j), we have the same equality on (—1,0) x {0}. 
Thus 

/ *? f)ll / _\ 

— dx i = log |g| - log (l + V 1 - 9 2 ) 

and 

pv / £ dxi=_ iog |?|+iog ( x+ ’ 

where the principal value can be ignored for exactly one of these integrals because there is no 
singularity. It follows that (1301) holds. □ 

Proof of Lemma [5] Let r £ (0,1) and consider the integral 

I r = |Vit| 2 da;. 

JR 2 + \B r (0) 

In order to estimate 7 r , we first study the image of R+\B r (0) under F _1 and then perform the 
change of variables x = F(w) = — l/cosh(iu). Recall identity (l25l) . which implies that 


£ S : Rew < — arcosh — 1^ | C F 1 (R^_\R r -(0)) 


1 


C £ S: Rew < ^ arcosh ^4 + 1 

where S is the domain defined on page [13] Recalling the function u(w) = u(x) for x = F(w), we 
see that 


Ir = 


iF-i-iR 2 \B r (0)) 


|Vm| 2 dw £ 




Z arcosh ( 4 - 1) , Z arcosh (4 + 1 


Thus 


4 - 2 < e 2/r / 7r + e 2Ir/v <4+2. 


As 0 < e 2 +/ 7r < 1, this means that 


4 - 3 < e 2/r / 7r <4+2 


and 


7T / 4 \ 7T / 4 

5 log( ? -3)<7 r < I log(- 3 +2 


In particular, there exists a universal constant C i > 0 such that 


|Vu| dx — 7rlog - 
K=\B r (0) r 


<Cir 2 , r £ (0,1). 


( 33 ) 


For b £ (—1,1) and r £ (0,1 — |6|), we have 

|Vub| 2 dx = 

Jr 


lR 2 \B r (b) 


l\*_ b (B+(b)) 


|Vw| 2 dx. 
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Thus we now examine the set c f>-_ b (.B+(&)). Since is a Mobius transform, it maps the semicircle 
d + B r (b) to a semicircle, which contains the points 

$_&(&-r) = --- J < 0 and <h_ b (6 + r) = ---— > 0. 

1 — b z + br 1 — b z — br 

Since <l>_ b (6) = 0, we obtain 

-®r/(l-b 2 +|b|r)(0) C ®~b(B+(b)) C | b | r )(0). 

It then follows from (l33l) that 


2 C-, r 2 

’ 108 - + T l°g(l - - |i|r) - (1 _y,_ |t|r)i 

< I \\7ub\ 2 dx 

Js.2_\B r (b) 

2 Ci r 2 

< ^ log - + ^ log(l - b 2 + \b\r) + (1 _ b2 _ |b|r)2 ■ 

Hence (l3lT) holds. □ 


Proof of Proposition [?[ The formula follows directly from the definition of W ±, the definition (IT5T) 
of E* d , Proposition [5j and the last two lemmas. □ 


2.3 The rescaled tail profile 

Since u* d is the expected limit of the rescaled stray field potential, assuming that condition (|4]) will 
be preserved in the limit, we can derive an expected profile p* d for the rescaled first component of 
the tails of the Neel walls. 

Using the unique solution u of the problem (fl9l) (I22f) . we first define the logarithmically rescaled 
tail profiles 

p(xi):=-PvJ ^-(t,0)dt ^ log ^1 + yjl - xfj - log|xi|, Xi e (-1,1) \ {0}, (34) 

and 

Vb(xi) ■= M 0< ^-b(^i) , xi € (-1,1) \ {&}• (35) 

Then for = u o <h_ 6 , we find = —p' b on (—1,6) x {0} and (6,1) x {0}, using the 

conformality of and the fact that is real on (—1, b) x {0} and on ( b , 1) x {0}. If we define 


N 


P 


\,d= Xi S (-1,1)\ Oat}, 


n= 1 


then we also have 


du* 


dx 


— = —(/x* d )' on (—1,1) x {0} except at the singularities a k , k = 1,..., N. 


Note that 


Pa n (a k ) = Pa k {a n ) = log 


1 + \A ~ g(a n ,ak) 2 

g{a n ,a k ) 


, n k. 


We need to examine the behaviour of p^ d near the points a n , as this will be important for 
determining the energy of the tail-core interaction. 
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Proposition 10. Set 


A n = In log(2 - 2 a 2 n ) + 7 kPa k M 

k^n 

for n = 1,..., N. Then there exists a constant C = C(a ) > 0 such that 


Ta,d(. a n i r ) A n 


In log - 
r 


< Cr 


( 36 ) 


for any r £ (0, p{a)\ and any n = 1,... ,N. 

Proof. First we study fib again for a fixed b. Obviously, for c £ (—1,1) with b ^ c, this function is 
smooth at c, and therefore there exists a constant C\ > 0, depending only on b and c, such that 

\p b {c±r) - p b (c )| < Cir 

for any sufficiently small r > 0. We also have, by (1M1) and (l35l) . the formula 


hb(b± r) = log 1 + Wl - 


(1 — b 2 =F br) 2 
Thus there exists a constant C 2 = C 2 (b) such that 


p b (b ±r)~ log(2 - 2b 2 ) - log - 

r 


+ log(l - b 2 =F br) + log ■ 


< C 2 r 


for r £ (0,1 — |6|). If we sum up for b = ak, k = 1,... ,N, then the conclusion follows. □ 

We will prove in Section [6721 that the function 


cos a 


Ha,d 


log 7 


gives an approximation for the profile of the first component m^ e for a minimiser m e of E e over 
M(a,d). In a standard Neel wall with rotation angle 2a (or 27r — 2a), the function mi >e would be 
required to make a transition from cos a to 1 (or —1) and back. If we superimpose several Neel 
walls, then this is no longer true. Instead, in an interval (a n — r,a n + r), up to a small error, ?ni j£ 
is now required, by (l36l) . to make a transition from 


cos a 


A n + J n log A 


to ±1 and back. We discount the term involving logA ; as it will be cancelled by a similar term 
elsewhere. The contribution of the remaining term to the energy (to leading order, rescaled by 
(logd) 2 ) is then —n 7 n X n near the point a n (see Remark [TUI below!. Thus in total, we have a 
correction term of the form 


N N / 

W 2 (a,d) = —7r^7„A„ = 7nlog(2-2a 2 ) + ^7„7 fc /x afc (a„) ] . 

n =1 n =1 \ k^n 

Using the above explicit expression for p, ak (a n ), we find 


N N 

W 2 (a,d) = -tv Y,^og(2-2a 2 n )-nJ2J2 Thin log 

n=l n=1k^n 


1 + y/l- pjak^n) 2 
ei.a k ,a n ) 
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This corresponds to the sum of the tail-core interactions described in Section [L3l Note that 


W 2 (a,d) = —2W\{a, d). 

We finally define W(a,d) := W\(a, d) + W 2 (a,d). That is, 

JV JV 

W(a,d) = -\Y, In log(2 - 2a 2 ) - - ^2 7fc7 n log 

n= 1 n= 1 k^n 


( 1 + a/1 ~ g{ak,a n ) 
y g{a kl a n ) 



3 The Euler-Lagrange equation 

Since we study the number inf M(a,d) E e for a given a £ An and d £ {±1}^, it is useful to consider 
minimisers of E e in M (a, d) (see Proposition [T]). The problem gives rise to an Euler-Lagrange 
equation, which is most easily expressed in terms of the continuous function ip : (—1,1) —> R with 
ip{— 1) = a and m = (cos ip, sin <p) £ TT 1,2 ((—1,1), S 1 ). Let u = U{m) be the function defined on 
page OH Then the equation is 

Du 

eip'\x i) = -—(a;i,0)sin(p(a;i) for X\ £ (-1, l)\{ai,..., oat}. (37) 

OX 1 

For the derivation of (1371) . let C £ Cq°(— 1,1) with £ (a*) = 0 for k = 1,... ,N be an infinitesimal 
variation of mi. Then 

(If |VC/(mi + tQ\ 2 dx ) = [ Vt/(?m)- VU{Qdx 
t=o \ 2 Jvt 2 + J 



Equation (1371) is now obtained with the usual computations. 

We use the abbreviation u' for the derivative of the trace of u with respect to (—1,1) x {0}. 
Then we have the following shorthand form of (l37l) : 

tp>" = i/siny? on (—1, l)\{ai,..., a^v}- 

We now analyse this equation. More precisely, we prove an interior l / E 2,2 -estimate for solutions of 
(EZD and we prove a Pohozaev identity. As a consequence, we eventually find that the exchange 
energy in the core of a Neel wall is of order 0(l/(log(5) 2 ). 

3.1 An interior W 2 ^-estimate 

We have the following interior W 2,2 -estimate for a solution of the Euler-Lagrange equation and for 
the corresponding stray field potential. This estimate will be used in Theorem 1221 below in order to 
find the specific profile of the Neel wall. 

Lemma 11. Let 0 < r < r' < R! < R. There exists a constant C > 0 (depending only on r' — r 
and R — R') such that the following holds true. Let e > 0 and set I = (-R, —r) U (r, R). Suppose 
that the functions u £ W 1,2 (Bj { ( 0) \ B r ( 0)) and (p £ W 1,2 (/) solve the system 

0 in B+(0)\B r (0), 

ip' sin ip on I x {0}, 

u sin ip in I. 


Au = 

du 

dx 2 


d_ 

dt 
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Further suppose that sin 99 ^ 0 in I. Then 


L 


I V 2 u| 2 dx + e 


(0)\B r ,(0) 


l(-R',-r')U(r',R') 


(■ <p"f + (v ')(1 + COt <£)) dxi 


< c 




dx i + 


S+(0)\B r (0) 


|Vu | 2 dx 


Proof. We first note that is smooth on I (see [14]) and therefore u is smooth in B^( 0) \ B r ( 0) up 
to the boundary I x {0}. Consider the function v = on B^( 0) \ B r ( 0). We have 


sin ip 


dv d 2 

and —— = — —(cos <p) on I x { 0 }. 

OX 2 


Let 77 £ C^°(S^(0) \ B r (0)). Since u is harmonic in B^(0) \ B r ( 0), integration by parts yields 

[ 77 2 |Vu| 2 dx= f f] 2 v —— q (cos(/?) dx\ — 2 [ r]vX7p-X7vdx 

Jb+{ o) Ji d x i J b+( o) 

= —e / 7y 2 -^—(<//'sin<p + (y/) 2 cost/j) dxi — 2 / rjvVr/ ■ Vv dx 

Jj sm(p J_b+( o) 

=—e f if(<p") 2 dx\ — e f 7 f ^ cot dxi — 2 f rjvVri-'Vvdx 

Ji Ji 3 ^s+(o) 

= — ej r) 2 {ip") 2 dx\ — ^ J t? 2 (^') 4 (1 + cot 2 y>) dxi + ^ J pr]'((p 1 ) 3 cot ip dx\ 


— 2 / r/uVr/ • Vw dx. 
ds+(o) 

By Young’s inequality, 

riri'(tp') 3 cot ipdxi < 


and 


Therefore, 


^ J d 2 <y) 4 cot 2 ip dxi + J (?/) 2 (</?') 2 dxi 

[ 77wV?7 • Vu dx < i / 7? 2 |Vd 2 dx+ f |V?7| 2 7; 2 dx. 

Jb+(o) 4 J B +(o ) J B +( o) 


lB+(0) 


'Bb(O) 


f 77 2 |Vu| 2 dx+e f rf 2 {{pi") 2 + (v?') 4 (l + cot 2 y>)) dxi 

Js+(0) Ji 

ijynv'f 


(38) 


< 12 / u 2 |V? 7| 2 dx + 4e 


dxi. 


As it is harmonic, we observe that that 


(9 2 u 

9x2 (9xi ’ 


so that 



dv 

dx\ 


dv 

8X2 


in B+( 0 )\B r ( 0 ). 


Choosing suitable cut-off functions 77 , we can now easily derive the desired inequality. 


□ 
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3.2 A Pohozaev identity 

As in the theory of Ginzburg-Landau vortices, a variant of the identity due to Pohozaev [28] is useful 
for our problem. For a function u : B+ (0) —> R, we use the notation 

d 0 u = t— r ■ V'tt. 

1*1 

Lemma 12. Let e > 0. Suppose that the functions u £ W 1,2 {B+ (0)) and p £ W 1,2 (— 1,1) solve the 
system 


A u = 0 
du 
dxi 


in B+ (0), 


= ip'simp on (— 1,1) x {0}, 


ep" = v! sin <p in (—1,0) and (0,1). 

Then for any r £ (0,1) , 

e f (p 1 ) 2 dxi = re(p'(r)) 2 + re(p'(—r)) 2 + r f (\X7u\ 2 -2(d p u) 2 \da. 

J-r Jd+Br(0) ' ' 


Proof. As u is harmonic, we calculate 

'1 


div -|Vn | 2 x — (x ■ Vu)Vu 1=0 in B{ (0). 


For any r, s £ (0,1) with s < r, it follows that 

0= [ 2 -r(d p uf) da- f (^|Vrt| 2 - s (d p u) 2 ) da 

Jd+Br .fOI 7 Jd+B.tO) 7 

f 


L 


s du du , 
Xi———— dx i 


dx\ dx2 

Using 00 and dun, we compute 
r s du du 


d+B s ( 0 ) 
du du 

Xi-R —a— dx 1. 

OX i OX 2 


Xi 


dx i dx 2 


dx i = e / X\p p dx i = 


Xi- — (p') 2 dx i 
dxi 


= y {T'(-r )) 2 - y (v?'(-s)) 2 - | J (p ) 2 dxu 


Similarly, 


l = ? (v ' (r))2 “ j M s » 2 ~ if. 


Hence 


(39) 

(40) 

(41) 


(p’) 2 dx 1 + ^ (</) 2 <f*i = ^ (^(v\r)) 2 +e(p'(-r)) 2 + ^ (|Vu| 2 - 2 (<9 p n) 2 ) da 

- | |d</0)) 2 + e(^'(-s)) 2 + ^ (|Vn| 2 - 2 (<9 p n) 2 ) do^j . 
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Since ip' G L 2 (—1,1) and |Vu| G L 2 (B + (0)), there exists a sequence s*, \ 0 such that 


s fc ( e(^'(s fe )) 2 + e(^'(-s fe )) 2 + [ |Vi 

V Jd+B Sk (0) 


I 2 da 


0 . 


Therefore, we obtain the desired identity. 


□ 


As a consequence, we obtain the following estimate, which implies that the exchange energy in 
the core of a Neel wall is of order 0(l/(log<5) 2 ). 

Lemma 13. Suppose that e G (0, |] and the functions u G W 1,2 (B+ (0)) and <p G W 1,2 (—1,1) solve 
the system ©D-HU). Let 


F = | f O ') 2 dx i + \ [ |Vtt| 2 dx. 
1 J -1 1 JBt( 0) 


Then 


r s 

/> 


dx i < 


Bf( 0 ) 
F 


lo S« 


l ■ 


If </>(0) G 7rZ, then 

Proof. It follows from Lemma IT2l that 


I s 

/ sin 2 pdxi < 85F. 
J-s 


(42) 

(43) 


1 1 


2 Js r J - 

Thus there exists an r G (<5,1] such that 


f (if') 2 dxi dr < F. 
J —r 


(p') 2 dx i < 


F 


Inequality (|42ll then follows immediately. 

For the second inequality, we note that for \x\\ < S, 


| sinyj(xi)| = 


rx i 

/ ip'(t) cos p(t) dt 
Jo 


1/2 


<[\xi\J o {v') 2 dt\ < 


2 5F 


elog 


1 — 


< 2 VF. 


(Here we use that fact that log 7 < 2 log | for 0 < e < 1.) Thus (l43l) follows immediately as well. □ 

Later we will prove estimates similar to (1421) and (l43l) even without making use of the Euler- 
Lagrange equation, but assuming suitable control of the energy instead (see Theorem[28]and Remark 

m- 


4 The core 

In this section, we study what happens near a single Neel wall, rescaled to unit size. Several technical 
difficulties arise in the analysis of the local behaviour of the magnetization, due to the nonlocal nature 
of the magnetostatic energy. For this reason, we first introduce a modified energy functional where 
the stray held is considered on a half ball (instead of R^_) with a Neumann boundary condition. 
We prove energy estimates (upper and lower bounds) for this modified functional together with 
statements about its behaviour under small perturbations of the boundary data. The analysis of 
the minimisers of this functional is essential because it yields good approximation results for the 
Neel wall profile (see Theorem [ 22 l) and for the core energy (see Theorem©]). 
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4.1 A modified functional 

For 7 £ (0, 2), we define the convex set 


m 7 = e W 1,2 (—1) 1) : m( 0) = 1, M±l) < 1 - 7}- 


We think of /z as the first component of the magnetisation near a Neel wall with transition angle 
2 arccos(l — 7 ). For /z £ M 7 , consider the convex set 

( 0 )) = {w £ W^(B+ ( 0 )): w(xi, 0 ) = p(x 1 ) for X\ £ (— 1 , 1 ) and w < 1 — 7 on 9 + _Bi(0)}. 

Clearly W*’ 2 (B+ (0)) yt 0 and there exists a unique function where the infimum 

inf f \S7w\ 2 dx 
wGW^ 2 (B+(0))Jb+(0) 


is attained, owing to the strict convexity of the Dirichlet functional. 
Now we define the following modified energy functional for p £ My 



^ ^ dx 1 + j- inf [ |Vtn | 2 dx. 
1 - h 2 «>ew4 2 (s+(o)) Jb+{o) 


Proposition 14. T7ie functional Ef admits a unique minimizer p £ M 7 , which satisfies /z(±l) = 
1 — 7 and 1 > /z > 1 — 7 in (—1,1). Moreover, if v £ (0)) is the unique function with 

f \Vv\ 2 dx = inf f |Vu >| 2 dx, 

Jb+( 0 ) weWfi 2 {B+( 0 )) JB+(0) 

then 1 — 7 < t> < 1 in B+(0) and v is the unique solution in W 1,2 (B+( 0)) of the boundary value 
problem 


0 

II 

<1 

in B+ ( 0 ), 

(44) 

M^i, 0 ) = p{xfi) 

for x 1 £ (- 1 , 1 ), 

(45) 

v = 1 — 7 

on 9 + i?i(0). 

(46) 


Proof. The direct method of the calculus of variations leads to existence of a minimizer /i £ My of 
Ef. Moreover, if /z £ My is a minimiser of Ef, then using a simple cut-off argument at 1— 7 and 1, we 
see that M±l) = 1 — 7 and 1 > p > 1 — 7 in (—1,1). This implies that the corresponding minimizer 
v of the Dirichlet energy satisfies 1 — 7 < v < 1 in B+ (0) and solves (l44l) (l46l) . Obviously, this 

boundary value problem has a unique solution in W 1,2 (Bt (0)). Since the function (aq, X 2 ) 1 —> , ; 

1 1 x 2 
for (aq,aq) £ M x (—1,1), is strictly convex, we deduce that Ef admits in fact a unique minimizer 

P £ My. □ 

Let pi £ M 7 be the minimizer of E] and v be the solution of (l44l) (l46l) . Since curlVM = 0, 
there exists a function u £ bF 1 , 2 (lM(0)) with V ± v = -Va. This is then a solution of 


A u = 0 in B+ ( 0 ), 

= -A on (- 1 , 1 ) x { 0 }, 
x • Vjj = 0 on dB+ (0). 


(47) 

(48) 

(49) 


Note that u is determined uniquely up to a constant by these conditions. On the other hand, 
given a function u £ W 1 ’ 2 (B^ ( 0 )) satisfying (1471) (141)1) . we can reconstruct a corresponding function 
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v £ 0 )) with 

the quantity 


—Vzt. It follows that the minimiser p of Ef automatically minimises 


e 

2 



( m ') 2 

i-a * 2 




|Vm | 2 dx, 


where u is determined (up to a constant) by (HT1) fill . Since p plays the role of the first component 
of the magnetisation near a Neel wall, then u roughly corresponds to a stray-field potential associated 
to n in ( 0 ). 

The Euler-Lagrange equation for the minimiser of Ef is therefore 


eip" = v! sin p in (— 1 , 0 ) and ( 0 , 1 ) 


(50) 


for any continuous function ip : (—1,1) —> R with p = cosy?, similarly to Section [3] 


4.2 Energy estimates 

We first prove the following preliminary estimate with a direct construction. This result is similar 
to {?] (see also d). 

Proposition 15. Let (3 £ (0,1). Then there exists a constant Co > 0 (depending on (3) such that 
for every 7 £ ((3,2 — (3) and every e £ (0, 5 ], 


inf e:<^L + . 


m; e 2log| (iogi) 2 ' 

Proof. Consider the function p : (—1,1) —► [1 — 7 ,1] given by 

^,) = l- 7 to8( r- + f' ) h 1 ° 8 f C[l- 7 ,1] for i, e (-1,1). 


log(l + S 2 ) — logci 2 


Then 1 + p > 2 — 7 and thus 


1 - (M^i)) 2 > 7(2-7) 


log(xf + <5 2 ) — log<5 2 
log(l + 1 5 2 ) — log <5 2 


Therefore, 


1 


a 


wr_ 

p 2 


■ dx 1 < 


2 7 2 


x\ 


7(2 - 7 ) log + 1 (*1 + S 2 ) 2 log (f* + l) 

2 7 2 +2 


dx 1 


S 7 (2 - 7 ) log y/^+1 ^ + 1)2 l0g( * 2 + 


dt. 


Define w : B+ (0) —> K. by 


for 0 < r < 1 and 0 < 9 < ir. Then 


w(r cos 9, r sin 9) = p(r ) 


[ \Vw\ 2 dx = - nl „ I \ 

(log^^Ti) " {T +S ) 


dr 


7T7 


2 (log V / JTT) 


L 


l+5‘ + _ p2 
5 2 t 2 


t — 6 2 , 7T7" 

dt < ’ 


!°g \ -ki + 1 
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Hence 


e:(p) < 


7T7 


2 log J 


3* + 1 


+ 


7(2 - 7 ) log i log V /J7I & + X ) 2 l0 S( i2 + !) 
which implies the statement of the proposition. 


dt , 


□ 


We can match the leading order term in this inequality with an estimate from below. Moreover, 
we obtain more information about the behaviour of the minimiser /i of E\]: in particular, we have a 
uniform W 1,2 -estimate for the difference between the rescaled stray field potential it log | (associated 
to /i via (S3 (ED) and the map 


(27,2:2) H > 7 ( arctan 


TTXl 

2|aq | 


(27,2:2) G 


First, however, we need some information on the regularity of solutions of the Euler-Lagrange 
equation, especially at the boundary. 

Lemma 16. Let 7 G (0,2). Suppose that ip G W l,2 {— 1,1) with cosip(— 1) = cos^>(l) = 1 — 7. 
Further suppose that u G W 1 ' 2 (B^ (0)) is a function such that CTl) (15(J1) are satisfied for p = cos<p. 
Then V« is continuous in B+( 0) and has a continuous extension to dB+ (0)\{0|. 

Proof. Let p = cos — 1 + 7 and consider the extension of u to R 2 . and of p to R given by 

u(x) = u 


and 


for x G R^_ with |®| > 1 


p(x 1) = —p — I for 27 G (—00, — 1 ) U ( 1 , 00). 


Then we have A u = 0 in 


and ^ = -p! on 


x { 0 }. Since p! G L 


loc \ 


we conclude 


that u'(-,0) G Lf oc (R) for any p G [1,2) with standard regularity theory. Thus by ED, we have 
tp" G L p (— 1,1), and it follows that p! is Holder continuous on [—1, —r] U [r, 1] for every r G (0,1). 
The extension is then locally Holder continuous in R\{0}, and using standard regularity theory once 
more, we conclude that Vw is continuous away from 0. □ 

Theorem 17. For any (3 G (0,1), there exists a constant C > 0 (depending on (3) such that the 
following holds true for every e G (0, |]. Suppose that 7 G (/?, 2 — /3) and p G M 1 is the minimiser 
of E, 7. Then 


Ef(p) - 


7T7 


2 log i 


< 


c 


( lQ g jY 


Let u G W 1,2 (B(( ( 0 )) be the solution of (l 47 l) (l 49 l) . Then 




b+{0)\b 5 (0) 


Vit(a;) — 


72; 


and 


Id+B 1 ( 0 ) 


x\ 2 log 
|Vu| 2 da < 


dx 


C 


L 


Bs( 0) 


|Vu| 2 dx < 


C 


( lQ g 


( lQ g 7) 


2 ' 


(51) 


(52) 


( 53 ) 
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Proof. We already know, by Proposition [l5l that there exists a constant Co > 0, depending only on 
(3, such that 


E7(v) < 


21og| + 


Co 

(i°gi) 2 


(54) 


We want to prove an estimate from below of the same order. To this end, we choose some <p £ 
W 1,2 (—1,1) with /i = cos ip. Consider the function £ : B+{ 0) —> R, defined in polar coordinates by 


£(rcos0,rsin0) 


7(0 - tt/2) 

i°gy 

yr (0 — 7 r/ 2 ) 

<5 logy 


if r > 6. 

if 0 < r < S. 


According to Proposition [HI we have 



= an.oMu)*! - J_' (2j^r + «*>■ 0)) 

Using Lemma fT3l and (f54l) . we find a constant C\ = Ci(/3) such that 


n'(xi) dx i. 


/ 77ra;i 


+ f(xi, 0) ) n'{x\) dx i 


l-s \2|xi| log y 
Moreover, by (TX71) (HU1) . we have 

J £(xi,0)n'(xi)dxi = 

Combining these estimates, we obtain 


< 


77T 


2 log ■ 


r s 

— I \<p'\\ sin<p| dx i < 

A JS 


Cr 


(log y) 


2 ’ 


[ f— — dx i= f V£ ■Vudx. 

/( —1,1) X {0} C'2-2 ./r+c 


/B+(0) 


7 2 7r Ci 


We compute 


/ V£ • Vu dx> , 

Ib+( 0 ) logy (logy)' 




and 


This implies 


|V£| 2 = 


( lo gy)‘ 


if r > S. 


[ |V£| 

Jb+< 0) 


9 , 7 2 7T 7 2 7T /I 7T 2 \ 7 2 7T 7 2 7T 

2 dx = -E-r + ---* I - + — ] < ' ' ' 


lB+(0) 

Hence, by (l54l) and ([55)l . 


logy (log|) 2 \2 24; logi (log±) 2 ' 


/ \X7u-X7f\ 2 dx <2Ef(n) 

JB+( o) 


7 2 7r | 2Ci + 7T7 2 Ci 


log 7 (logy) 2 (log 


(55) 
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where C 2 = 2Cq + 2C\ + 7T7 2 . It follows in particular that 


/ |Vrt| 2 dx < 2 f \X7u-X7£\ 2 dx + 2 [ |V£| 2 cfa< 

Jb+( O) Jb+(0) JB+(0) 


iBfiO) 


and 


L 


S+(0)\S 5 (0) 


\7u(x) — 


"IX 


\ x \ 2 log | 


BJ( 0) 


dx < 


2 C 2 + 2 tt 7 2 

( lo s |) 2 


C 2 


(logi)- 


since V£ = rulr—r if r > <5. Furthermore, we have 
M 2 log i 


[ \X7u\ 2 dx= [ |Vit-V£| 2 dx- f |V£| 2 dx + 2/ X7f-X7udx 

Jb?(o) Jb+(o) Jb+( 0 ) Jb+io) 


|Vrt| 2 dx = 

'B + ( 0) JB+( 0) 

7 2 7t 7 2 7t + 2C\ 

“log| (logi) 2 ■ 

If we combine this with (1541) . then we obtain 

rl (m ') 2 


>B+( 0) 


: [ d&i < 2Ef ht) - [ |Vit| 2 dx < 

J-l 1 — P «/.B+(0) 


2 Co + 2Ci + 7 2 7r 


Thus we have proved inequalities (15T1) and (15^1) . 

Finally, we apply Lemma IT21 and let r 1. Then by Lemma fl6l we obtain (15^1) . 


□ 


4.3 Behaviour under small perturbations 

Next we want to understand how the number infjw Ef changes when we vary 7, and how the 
energy changes when we perturb the boundary condition (1491) . In particular, we prove the following 
statements. 

Proposition 18. Let fd £ (0,1). There exists a constant C(/3) > 0 such that for all 71,72 £ (/?, 2-/3) 
and every e £ (0, 5], 

iof E? - EJ‘ - JUL + 

M -V2 2 logi M 71 2 logi (log i) 2 

In other words, the function 

2 

g : (0,1) -4 R, 74 inf Ef - l t (56) 

m 7 2 log | 

is locally Lipschitz continuous with Lipschitz constant of order 0(l/(log<5) 2 ). 

Proof. First note that it suffices to prove the inequality when 172 — 711 is small, for otherwise it 
follows from Theorem [171 

Let hi £ M 7l be the minimiser of Ef 1 and let v\ £ W 1,2 (B2 (0)) be the corresponding solution 
of (El) (El) in Proposition [TT] Define 

M2 = —(Ml - 1) + 1 
71 

and 

v 2 = — {vi — 1) + 1. 

7i 
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Then G M 72 and u 2 is the unique solution of (l44l) (l46l) associated to M2, so that, by PronositionUdl 




dx. 


We have 


and 


Hence 


1 + M2 


— (1 + Mi) + 
7i 


2(7i ~ 72) 
7i 


1 - M2 = —(1 - Mi)- 
7i 


i-Mi 


(m'i ) 2 

1 -M 1 

1 I 2 ( 71 - 72 ) ■ 
72(i+ah) 


Under the assumptions of the proposition, we have /3 < 71, q 2 < 2 — /3 and /3 < 1 +mi < 2 throughout 
(—1,1). Therefore, we have a constant C\ = C\((3) such that 


1 


1 , 2 ( 71 - 72 ) 

72(l+/2l) 


< 


72 


-7? 


■ Ci|7 2 


7i|> 


provided that |q 2 — 711 is small enough. 

We now have 

£?(M 2 )< (4+Cil72-7i ^ }^dx 1 + £ f |Vr>i| 2 dx. 

V7i / 2 1 — Mi 2 7 f (0) 

Combining this with (1FT1) and (l52l) . we obtain the inequality 


5(72) < 5(71) + ^2 a g(7i) + 


C2I71 - 72 1 

( lo §i) 2 


for a constant C 2 = C 2 (/3) by Theorem[l71 where g is the function defined in (l56l) . Finally, we know 


that 


3(7i) < 


C 3 

( lo s |) 2 


for another constant C 3 depending only on ft (by Theorem [T7] again). Hence we obtain the desired 
inequality. □ 


Remark 19. Recall the discussion in Section \2.3t Consider a Neel wall at the point a n such that 
nii makes a transition from 1 — 7 n to 1 and back. If we modify the wall, changing j n to 7„ — £ n 
with ( n = then the change of the energy (to leading order, rescaled by (log<5) 2 ,) is 



( inf 

V M 7„-C„ 


- inf Ef n 


-7T7nA„ + 0(1) as e \ 0. 


This is the phenomenon that leads to the core-tail interaction term in the renormalised energy. 
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Lemma 20 . Let Co > 0, (3 £ (0,1) and q > 2. Then there exists a constant C > 0 (depending only 
on Cq, (3 and q) such that for any 7 € (/3, 2 — (3), any e £ (0, |], and any 77 £ (0, Co), the following 
holds true. Suppose that p £ M 1 and u £ W 1,2 (B(~( 0)) with Au = 0 in B+( 0) and = — 7 / on 
(—1,1) x {0}. Suppose further that 


and 


Then 


||z • Vu|| L <i(a+ Bl (o)) < 


Vu — 


7 X 


x \ 2 log I 

e I 1 (MO 2 


I|Vm|| L 2 (b + (0)) < 

+ fn\\ n.fnw AW o s 


2 J-x 1 - T 2 


dx 1 + 


L*(B+(0)\B,( 0)) 

1 r 


2 Jb+{ 0) 


|Vu| 2 dx > inf Ef — 


Cv 




(i°g 7 ) 


2 ’ 


(57) 


Proof. Consider the function w £ W 1,2 (B(~(0)) with Vw = X7 ± u and u>(xi,0) = p(x±) for x\ £ 
(—1,1). Let Sg = {(cost, sin t): t £ (0, 9)} for 0 < 9 < z. Note that 


f 9^ 7 r 

(cos 9, sin 9) — p(l) = / x-Vuda< - 1 —, 9 £ (0, —). 

JSe log i 2 


Similarly, we find that 


Q t 7j- 

u)(— cos9, sind) — p(— 1) < - - —, 0e(O, —). 

logy 2 


Thus if we define 


and 


9{x) 


log 7 


w = w — g, 


then we have that w < max{/r(±l)} < 1 — 7 on dB(f( 0) (because 2sind > 9 for 9 £ (0,-|)). 
Moreover, w £ W X,2 {B^( 0)). Indeed, we have 



L 


B\ ( 0 ) 


|Vu| 2 dx 


'L 


B\ (0) 


Vit • g dx + 



|V.g| 2 dx. 


For any p £ [1, q), we have a constant Ci = C\(p) such that 

CiT? 


l|V.g|| 


< log i- 


For p > 2, inequality m gives another constant C 2 = C 2 (p, Co) such that 

C 2 


I|V'u|Ilp/(p-i)(B+(0)) < 1 ' 

We conclude, using Holder’s inequality, that 


L 


B+(0) 


„ . C\C 2 7J 

vu • V gdx < 


(log 7)' 
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using some fixed p £ (2, q). Therefore, 


f \Vw\ 2 dx< f \S7u\ 2 dx-\ - 

Jb+( o) Jb+( 01 floeri 


B+(0) 


C 3 i! 

(log l) 2 


for a number C 3 > 0 that depends only on q and Co- Now the statement of the lemma follows. □ 

Corollary 21 . Let (3 £ (0, |), Co > 0, and q > 2. Then there exists a constant C > 0 such that 
for any 7 £ (2/8, 2 — 2/3), any e £ ( 0 , ^], any p £ ( 0 , Co), and any ( £ (—Co, Co), the following holds 
true. Suppose that p £ W 1,2 (—1,1) with p(±l) < 1 — "f +-r^-r and p( 0) = 1. Let u £ W 1,2 (Bt (0)) 

1o S s 

be a function with A u = 0 in B+ (0) and = —p! on (—1,1) x {0}. Suppose that 


\\x ■ Vu|| i 9 (a+ Bl (o)) < 


lo 8 i 


and 


Then 


Vu — 


72 ; 


\ x \ 2 log 7 


< 


Co 


II Vw IIl= ( b+ ( o)) ^ j 1 • 


L 2 (B+(0)\B s (0)) 


nl id') 


!\2 


2 7-1 1 - C 


■ dx 1 


\j + |V«| 

1 Jb+( 0 ) 


TT"/( + Cp 


c 


2 dx > inf Ef - 7 -=■• 

M -v (log ±) (log i) 


Proof. Set 7 = (J log g. We have p € M 7 _ 7 , thus for e small enough, we can first apply Lemma [21)1 
for 7 — 7 £ ((3,2 — (3) instead of 7 . We obtain 


nl (/*') 


/\2 


<ixi 


f \Vu\ 2 dx> inf S , 7 " 7 - 

B+(0) M r-S (log±) 


2 7-i 1 — d 

for a constant C\ = Ci(/3,Co,q). From Proposition IT51 it follows that 

7T7C C 2 


inf El -1 — inf Ef > — - 
m 7 _, m 7 (log I) (logy)' 


for another constant C 2 = C 2 (/3, Cq). Now it suffices to combine these estimates. 


□ 


4.4 The profile near the core 

Minimisers of Ef give a good approximation of the behaviour of Neel walls near the core. When we 
let e \ 0 , then, after renormalisation, we have convergence of those minimisers to a specific profile. 
More precisely, we have the following. 

Theorem 22. For e £ (0, and 7 £ (0,2), let p t £ M 1 be the minimiser of E]. Let u e be the 
corresponding solution of m AMD with 


L 


: dx = 0 . 


B+( 0 )\B 1 / 2 ( 0 ) 


(58) 


Suppose that 0 < r < R < 1. Then, as e \ 0, 


u e (x) log i —*• 7 ^arctan ^ — 2\ 1 \ ) wea ^V W 2 ’ 2 (B^(0)\B r 


m 


and 


(p e (xx) - I + 7 )log- 


7 log T —r strongly in W 1 ' 2 ((—R, —r) U (r, R)). 

Fil 
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Remark 23. For a related problem (concerning Neel walls in the presence of an anisotropy but 
without the confinement to an interval), Melcher \2f\ (25) proved similar results on the profile of 
Neel walls (with methods different from ours). 

The proof of Theorem [22] relies on Lemma [TT] Therefore, we first need to show that the as¬ 
sumption sinc^ e 7 ^ 0 in that lemma is satisfied for a function <p e with p e = cos ip £ , at least away 
from X\ = 0, for sufficiently small values of e. Because we will use this result in a slightly different 
context as well, we formulate it more generally. 

Lemma 24. There exists a constant C > 0 such that the following holds true. Suppose that Cl C 
is a Lipschitz domain with (—1,1) x {0} C dCl and let v be the outer normal vector on dCl. Let 
e £ (0, |]. Suppose that p £ W 1,2 (— 1,1) and let u £ W 1,2 (Cl) be a solution of 


A u = 0 in Cl, 


— = -p! on (- 1 , 1 ) x { 0 }, 
UX2 


v • Vu = 0 on dCl\{{— 1,1) x {0}). 

Let x\ £ (<5 — 1,1 — (5) and set 

S = (-B/(a;i,0) nfi)u{j/eUn B%( 0 ): y 2 > \yi - El I or \y\ £ (|, 2 )} . 


Then 



1/2 


Proof. Define 


iff) 


0 


1 

X\—t 

28 


if t < x\ — 8, 

+ 5 if X\ — S < t < x\ + S, 
if t > x\ + S. 


Furthermore, let 



and define x '■ R+ —> K by 


x(xi + r cos0,r sin#) = (1 — rj(6))f(xi + r) + rj(9)tf)(xi — r). 


Finally, let 



Then Cl n supp Vy C E and 



for some universal constant C\. 
We have 
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Moreover, 


/ I rxi+S 

x{t, 0)/j,'(t) dt < / |/z'| dt 

-1 J x-\—8 


rx i+5 


1/2 


< 26 


' x\—6 


(//) 2 dx 


Using the boundary value problem for u , we find that 

-l 


J x(t,0 )n'(t)dt = Vx-Vudx < ^Ci log i J \Vu\ 2 dx 


1/2 


Combining these estimates, we obtain the desired inequality. □ 

Proof of Theorem\EE We choose ip e : (—1,1) —> [0,7r) such that </> £ (0) = 0 and = cosip e . Using 
Theorem El and Lemma l24l (applied for f 1 = B+( 0)), we see that for any r, i? £ (0,1) and any 
sufficiently small /3 > 0, we have |/x e — 1 + y| < /3 and | sin</? e | > /3 in (— R, —r ) U (r, R) if e > 0 is 
small enough. Therefore, we can apply Lemma [11] in (— R, — r) and in ( r,R ). 

By Theorem [T71 the Poincare inequality, and (1581) . the functions u e log | are uniformly bounded 
in the space W 1 ’ 2 (i?( , “(0)\i3 r (0)) for all r > 0. Hence, by Lemma fill they are uniformly bounded 
in W 2 ' 2 (B^(0)\B r (0)) whenever 0 < r < R < 1. Hence there exists a sequence e*, —>■ 0 as k —> oo 
such that we have weak convergence 

Ue k log t ^ w 

Ok 

in W 2 ' 2 (B~l i (0)\U r (0)) for all r,R £ (0,1), where 6k = Cfc log The limit 


w€ p| W 2 ’ 2 (H+(0)VB r (0)) 

0<r<R<l 


is harmonic in £?+(0) and satisfies the boundary condition x ■ Vic = 0 on d + Bi(fS). According to 
Lemma E we also have 


lim sup 

k—> oo 



< OO. 


It follows from (l50ll that 


lim sup 

k—t oo 




(u tk ) 2 dxi 


< oo. 


Thus w is constant on (0,1) x {0}, and we can prove the same on (—1,0) x {0}. 
Define 

, xeB+(0). 


Since 


L 


B?( 0) 


Vic — 


X2 \ 

TTXl \ 

Xi) 

2\xi\)’ 

_L 

2 

M 2 

dx < oo 


by Theorem fl7l it follows that u> £ W 1,2 (B^ (0)). We conclude that u>(-,0) £ H 1//2 (— 1,1). More¬ 
over, the trace ic(-,0) is constant on (—1,0) and on (0,1). But iJ 1 / 2 (— 1,1) does not allow any 
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jumps; hence w{- , 0) is in fact constant on (— 1 ,1). We also have Aw = 0 in B+ (0) and x ■ X7w = 0 
on (? + i?i(0). Thus it follows that w is constant in B+( 0). Because of (l58|) . we have w = 0. That is, 


<x) = 7 (arctan (|) - ^|j) 


B,+ (0). 


As the limit is thus independent of the sequence £&, this implies the first claim. 

We have 

i 1 du *r n\ v 7 

log 7 7 T ~( Xl ’° ^ — 

0 OX 2 Xi 

strongly in L 2 (—R,—r) and in L 2 (r : R). But since /x' = —it follows that (/z e — 1 + 7 ) log y 
converges strongly in W 1,2 (— R,— r) to A_ — 'y log |xci| and in W 1 , 2 (r,i?) and to A + — 'y log |cci| for 
two constants A_,A+ G [— 00 , 00 ]. It remains to determine these constants. 

Choose a function x € Cq°(Bi(— 1 ,0)) with (—1,0) ^ suppVx- Then 


lim log 

e\o \ 


U 


B+(0) 


Vx • Vw e dx\ = 7 




11 X{xi) 

-1 X\ 


dx 1 . 


On the other hand, 


/ Vx • Vm £ dx= n' t xdx 1 = - (jif - 1 + 7)x' etei 

Jb+(o) J -1 J-i 


by an integration by parts. Hence 


lim 

e\,0 


log \J Vx-V« e da:J = J (7 log |a:i| - A_)x'otei 


— x(— 1 ) 0 )A_ — 7 J 

and we obtain A_ = 0. Similarly we show that A+ = 0. 


X(xi) 


_i Xl 


dx 


1, 


n 


4.5 The core energy 

We can now determine the values of the function e in Theorem O albeit not explicitly. They arise 
as the limits in the following result for 7 ± = 1 =F cos a. 

Theorem 25. For any 7 G (0,2), the limit 

106 7) 

exists. 

Definition 26. The function e : {±1} —> K. is defined by 

e( 1 ) — Ci+cosa and e(l) — ei_ cosa . 

Proof of Theorem 1,2,51 Define 

^>=H) 2 s 
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Fix e > 0 small enough and choose a number R £ (1, |). Let fi £ M 1 be the minimiser of EJ and 
let v £ W 1 ' 2 ^^(Q)) be the solution of (TH1) (11()1) . Define 


and 


fTlog^ 


l*i| 


A(*i) = < 


log? 

logi? 
log 


+ 1^7 


if ji < |*i | < 1 , 


1 - 


y) i) + t^t if |*i| < jz, 

\) !°g I 


v[x) = < 


VogT^y 
l°g A 


if ^ < |*| < 1 , 


+ ^ ifi*i<i 

V log \) log I 


Then we have jl £ M 1 and v £ W~ :2 (B + (0)). Moreover, by Proposition fl4l 

1 + fi(x i) > 1 + fi(Rxi) > 2 — 7 > 0 


and 


1 - p,{x i) = ^1 - (f “ M-R*i)) > o 


for < -h. Therefore, we compute 


i/« {ji , f 


L f 2 

R J-l/R 1 


-1/R t M 


■ dx i <11 — 


logi? 

log? 


1 (m 0 2 

-i 1 - M 2 


diri. 


Furthermore, 


K/r( 0 ) 


|V-D | 2 dx = ( 1 — 


logi? 

log J 




|Vu | 2 dx. 


We also observe that for x\ £ [—1, —-jj) U (j=r, 1], 


1 - fL 2 = 7 1 - 


lo %\k ' 

log 7 


2 — 7 + 


7 lo g^\ 


log 




Hence 


(AO 2 


■ dx i < 


2 7 


(-1 ,-i/K)u(i/fl,i) 1 - (i 2 (2 - 7) log y Ji/R x 2 (log i - log 


r f 


dx 1 


27 


(2 - 7 ) log \ 


l logi 


ds 


A - logs 


Define 


Finally, we compute 


g(R) = 


ilf: 


ds 


(2 - 7 ) log g J 1 log t - log s 


L 


2 , _ 7T7 ^g R 


|Vii| 2 dx = 

b+(o)\b 1/r (o) (log±) 
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Let 


Then we have 


i=~e(R) = -( 1- 


E 1 (A) < 


f \2 


g f 1 W) 

2 y_! i— 

log R 

toil 


■ daq 


2 ifl+(0) 


log-R\ 

log * y 

|Vw | 2 dx 


< e. 


< 1- 


ei m 


tg{R) f 1 _ logR\ | 7T7 2 log R 


R 


+ 


iog|; 2 (log *)" 


It follows that 


/e ) <( iog I) 2 (( i -|||) 2 £ 7 ( , ) + 


: g(R) 

R 


1 - 


k)gl?\ 7T7 2 logi?\ 7T7 2 1 


where 5 = e log i. Since we have equality for R = 1, we can use this inequality to estimate the 
left-hand superdifferential 

/i(e) = liminf «£WH. 

s/'t s — e 


Indeed, note first that 


and 


de 

dR 


= -e 1 + 


l 


1 


log* 


d6_ 

dR 


R= 1 


= e 1 - log 


1 + 


log* 


The above inequality therefore implies that for all e G (0,e 2 ], 




b6 - ii(i+ i^ 


- iog - ) E2 (/r) 




= (loglog - - 1 


1 -A ( 2E 7M 


V log* 2 log* log*/ 2-7 


1 

log 7 

67 


< 


C log log * 

(log*) 2 


for a constant C = C{ 7 ) by Proposition [15] Hence 

Clog log* 




: (l°g|)' 


Note that 


Thus if we denote 


f 


log log* 


6 (log *)" 


“77 dc < 00. 


e 7 = liminf /(e), 


e\,0 
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then for any > 0 we can find a number eo > 0 such that 

/( e o) < e 7 + ^ 

and at the same time, 

for any e £ (0, eo]. It then follows that 

/(e) < /(eo) ~ J f- (s) ds<e 1 + p. 

Hence we have in fact 

e 7 = 1 i n i/( e )’ 

as required. □ 


5 Several walls 

We now consider a given a £ Hjv and d £ {±1}^ and we study magnetisations m £ M(a,d). In 
particular, we want to estimate inf M ( a .d) E e and derive some inequalities for the magnetisation and 
the stray field in terms of the energy excess 

EJm ) — inf E e . 

M(a,d) 


5.1 Upper bound for the minimal energy 

The purpose of this section is to prove the following upper bound for the energy by a direct con¬ 
struction. It is a generalization of Proposition [15] to configurations with several walls a £ An- 

Proposition 27. For every R £ (0,1] there exists a constant Co > 0 such that for all e £ (0,1], 
a £ An with p{a ) > R, and d £ (±1} W , the inequality 


holds true. 


inf E e < 
M{a,d) 


2 M 


N 

^2(d, 


n= 1 


cos a) 2 


Co 

(log *) 2 


Proof. We will in fact prove a more explicit estimate, showing that there exists an m £ M (a, d) 
with 


EliK-cosa ) 2 


E e (m) < 


2 log 


R 2 


7r + 


t~ 


1 


sin 2 a log * Uoo {t 2 + l) 2 log(* 2 + 1) 


dt 


(59) 


This will clearly imply the statement of the proposition. The proof is similar to the proof of 
Proposition [15] Let 7 n = d n — cos a. Define 


/(2b) = 


log(:r 2 + (5 2 ) — log(l? 2 + S 2 ) 
logd 2 — log(i? 2 + S 2 ) 


and mi : (—1,1) [—1,1], given by 

mi{xi) = 


cos a + "/ n f(xi — a„) if x\ £ ( a n — R,a n + R) for n = 1,..., N, 
cos a else. 
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Then there exists a function m2 '■ (—1,1) —> [—1,1] such that to = (mi,m2) £ M(a,d). Suppose 
that X\ £ (a n — R,a n + R). If d n = —1, then 1 — mi(xi) > 1 — cos a and 

1 + mi (xi ) = (1 + cos a) (1 - f(x 1 - a n )) . 

If d n = 1, then 1 + toi(xi) > 1 + cos a and 

1 - mi(x\) = (1 - cos a) (1 - f(x 1 - a n )). 


In both cases, 

1 - (mi(ii)) 2 > sin 2 a (1 - f(xx - a n )). 
Thus as in the proof of Proposition [15] we can estimate 


L 


CL n +R 


27 1 


r-R 


xf 


| m | dx 1 < — - , / 2 \ 

a «~ R sin 2 a log yj ^ + 1 *~ R ( x l + l $ 2 ) 2 log (fj- + lj 


dx 1 


< 


27 ^ 


t z 


Summing over n, we find 
r 1 


\m' I 2 dt < 


5 sin 2 a log y/£ + 1 ( 42 + !) 2 lo s(* 2 + X ) 

2Eli7 2 r ^ 


dt. 


6 sin 2 a log ^f + 1 ^ + 1)2 log ^ + ^ 


dt. 


Now consider the function u = U(m) as defined on page [3] Since curlV^u = 0, there exists a 
function v £ IT 1 ’ 2 (R^_) such that Vv = \/ ± u in R+. Since this means that v' = m' x on (—1,1) x {0}, 
we can choose v such that v = m\— cos a on (—1,1) x {0}. Then we also have v = 0 on (— 00 , — l)x{0} 
and on (l,oo) x {0}, and of course Av = 0 in R+. Furthermore, the function has finite Dirichlet 
energy, and it follows that it is the unique minimiser of the Dirichlet energy under these boundary 
conditions. 

Define w : —» M by 

w(a n +rcosd, rsin8) = m\(a n + r) — coscn for 0 < r < R and 0 < 6 < n, n = 1,... ,N, 


while w = 0 in flp(a). Then we compute, similarly to the proof of Proposition 1 151 that 


L 


B+(a n ,0) 


|Vu>| 2 dx = 


7r 7r, 


log 


R 2 


+ 1 


R r 3 


0 (r 2 + (5 2 ) 2 


dr 


7r 7n 


2(log't/f+ 1 


/ 


R+ ‘^ dr< _ U 


lo S \/^ + ! 


Hence 


In particular 




|Vu;| 2 dx < 


E N 2 

n=l T 


log 


£ + 1 


[ \\/u\ 2 dx = [ |Vn| 2 dx < [ \Vw\ 2 dx < 

J R 2 + J j R’ logy/^ + 1 


If we combine these inequalities, then we obtain (1591) . 


□ 
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5.2 Stray field estimates 

The following is the main result of this section and one of the key ingredients for the proof of 
Theorem 0 We recall the function u* d from Section [2j solving A u* d = 0 in and ^ = 0 on 
(—oo,—l) x { 0 } and ( 1 , oo) x { 0 }, and with a piecewise constant trace on (— 1 , 1 ) x { 0 } given by the 
values 

( N n 

y, (dk - cos a) y(d k - cos a) 

k=n -\-1 k=1 

It has the property that for n = 1,..., TV, the function 



X ^ K,d( X ) ~ ( d n 


cos a) 


^arctan 


X2 

X\ dn 


7 r(xi - a n ) 
2\x\ - a n | 


is harmonic in LTt , ( a ni 0) and constant on ( a n — p(a), a n + p{a)) x {0}. Standard elliptic estimates 
then imply that this function is smooth near (a,j,0). In view of the energy estimates (1301) and (1311) . 
we can make more quantitative statements as well: if p(a) > R > 0, then 


du a,d , . , (d n ~ COS a)x 2 

-(a:) + --- 9 

OL Cl (xi - a n y + X2 


du* 


dx 2 


(dn - cos a)(xi - a n ) 
(xi - a n ) 2 + x\ 


< C 


(60) 


for x € Bp( a ) ( a n , 0), where C = C(a, N , R). 

Theorem 28. For any R £ (0, and Cq > 0, there exists a constant C\ > 0 such the following 
holds true. Let a £ An with p{a) > R and d £ {±1}^. Set 


N 

T = y^(d n — cos a) 2 . 

71=1 


Suppose that e £ (0, |] with S < R and m £ M(a, d) with 


E e (rn) < 


7tT 


Co 


21 °gj (log I)" 


Let u = U ( m ) be the function defined on page 0 Then 


(61) 


and 



(62) 


(63) 


This statement is somewhat similar to Theorem ll7l The main difference, apart from the fact that 
we consider several Neel walls here, is that we only assume a suitable bound for the energy, whereas 
in Theorem EH we study minimisers. Before we prove the theorem, we establish the following 
auxiliary result. 


Lemma 29. Let s > 0 and p £ W 1:2 (—s,s). If p{ 0) = 1 and |/u| < 1, then 



dx i < 2 s 



(M') 2 

1 — /i 2 


dx i- 
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Proof. By the Cauchy-Schwarz inequality, we have 

,/\2 \ 1/2 


Since /r(0) = 1, any x\ £ (— s,s) will satisfy 

1 - (m(zi)) 2 = -2 [ dt 

Jo 


1/2 


(64) 


< 2 


r x i !/ 2 /-xi 

/ / M 2 (i- M 2 )dt 

J o 1 — M Jo 


1/2 


Integrating over (—s, s), recalling that |/k| < 1, and using the Cauchy-Schwarz inequality, we obtain 

/ s / r s r Xl (a') 2 \ 1 ^ 2 / \ 1 ^ 2 

(1 — n 2 ) dx\ < 2 ( J J -- 2 dt dxij ( j J (1 — [i 2 )dt dx\\ 


(£aUru (i -^r 


< 2s 


which leads to 


f (1 — /i 2 ) dx\ < 4s 2 f 

J —S J —i 


/\2 


(/*') 

l-s 1 - /l 2 

Combining this with (l64l) . we obtain the desired inequality. 


dx\. 


□ 


Proof of Theorem\2f^ It is clear that it suffices to prove the inequalities for small values of e. We 
modify the functions u* d as follows: for a fixed s £ (0,1?], let € W b 2 (R2_) be such that 


60*0 = 


K,d 0*0 


for x £ n s (a) and 

6 (fl n + r cos 0, r sin 6) = 


log I 

ru* d (a ra + s cos 0, s sin 6) 


s\og\ 

for 0<r<s, 0<6*<7 t, and n = 1,..., N. Then we have 


0-D 


(Jn—1 “l - &n 

2!ogi 


f . l0 7 tT log - + Cl 

/ |V6| 2 dx< 2 

JrI (log i) 

for a constant C\ = C±(a,N,R) > 0. This follows from (1501) . (13TT) . and (RiUl) . 
We observe that 

ttT f 1 KJx i,0 ) 

i- T= —t-1— m^x^dx i 

log j J_i log j 

/*i />: 


(65) 


= / £ s (in,0)TOi(a;i)otei - / 60*4,0)- 


;, d (^i,o) 

log I 


to( (aq ) da:i 


f V£ s • S7udx — f f 60 * 4 , 0 )- a '^ / ^ m'^aq) daq. 

•/Mi J-l V log A / 
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We have a constant C 2 = C 2 (a, N,R) > 0 such that 




'-I 


Thus by Lemma [ 


/:,( 


6(®i,o) - 


We conclude that 


7lT 


< 


1Q g{ 

<A X i’°) 

log* 

2 C 2 s 


m[(xi) dx 2 


n N 


oa n -\-s 


log 


jm^l dx 1. 


(5 n=l J a « 


mi(x 1) ds] 


< 


2C 2 s 


elog 


T ( 2 - B e(m) - ||Vu||| 2(R 2 


r (2 E e (m) - IIV-uH^^)) + / 2 V£ s • Vitcfa;. 


log | e log i 

Using the Cauchy-Schwarz inequality and (1651) . we obtain 


ttF 2C9S / \ i/7rr log i + Ci 

i„7is Tijji ( 2£ -(“> - l |v “lh(«j)) + — iogi —» v “ll^cy 


( 66 ) 


(67) 


We want to use this inequality to prove (l63l) first. For this purpose, we choose s £ (0,1?] such 


that 


II^ u IIl 2 (r 2 ) 


7rr 


2 Ci 


logi (logi) 2 ' 

This is possible whenever e is sufficiently small because of m- Then (Ifni) and (16T1) imply 


( 68 ) 


ttT 2C 2 s / ttT t rT 2C 0 2Ci \ 

logi " elogi ylogi logi + (logi) 2 + (logi) 2 / 

y^r 2 T 2 (log i) 2 - CiTrTlog i -~2Cff 
+ log i log i 

2 C 2 s f 7rT ttT 2C 0 2Ci \ 

- elogi l^logi logi + (logi) 2 + (logi) 2 / 

+ *L(l __ 

log* \ 2^Tlogi ^2 (logi) 2 ; 

because y/1 — f3 < 1 — y for B £ (0,1). That is, 

£l + c 'i , <ig» ( toe 1_ tog l + 7c,iogl + 2 c,bgi y 

2ttF 7r 2 T 2 log * e log * \ S s S 5 ttT log * ttT log * ) 

In particular, there exist certain constants C3,C4,Cs,C6, all of them positive and depending only 
on a, N, Co, and R, such that 


C3 + 


C 4 < £ 
log i d 


log - + 
s 


C 5 log i C 6 log y \ 

logi + logi J 


(69) 


We want to use this inequality to show that there exists a constant C7 = Cr(a, R , TV, Co) > 0 such 
that C7S > 6. Then (l68l) follows immediately from (1 68 1) . because the right hand side is increasing 
in s when log * > 4 Ci/( 7 tT), which is the case for e small due to (ICTTl) . 

To this end, let c = min{l, We distinguish three cases. 
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Case 1. If s > cS, then the claim is obvious. 


Case 2. If s < c6 and 

1 £3 X] 1 

slog - > ^-0log -j, 
s 2C 5 0 

then it follows that s > 5 2 , provided that e is small enough. (Otherwise we would have an 
immediate contradiction to the assumptions for this case.) Hence log i < 2 log j and 


s > 


C 3 S 
4 C 5 


Case 3. If s < c5 and 


then we obtain 


, * C 3 ,, 1 

slog - < 7^-0 log T: 

s 2C 5 0 


C 


a 


~ + A ( l0g T 


5 , C 6 log i 


log I 


from m- We also have log i > log | (since c < 1). Hence 

!M< c <A 

dlogi - - 4C 6 - 


Hence 


which implies the claim. 


C 3 s S 
~r < t log-, 
4 0 s 


This concludes the proof of (1551) . 

Now we go back to inequality (1661) and use it for s = S in order to prove (1621) . Since we now have 
©sd, the inequality implies that 

f VZ s -Vudx >-^-—^-2 

Ul log 7 (logi)- 

for a constant Cg = Cg(a, TV, Co, R). Hence 

/ |Vu- V^| 2 dx < 2E e (m) — 2 / V^-Vudx+f |V&| 2 ds 
Jr 2 + Jr 2 + Jr 2 + 


(70) 


< 


2Co + 2Cg + Ci 

( lo g|) 2 


by (15T1) . (1551) and (1751) . Since £5 coincides with u* d /log| in fig (a), this finally implies (1551) . □ 

Remark 30. Once we have (ED and ©3D in Theorem \28[ we can also derive the inequalities 

C 


z J (ip') 2 dx 1 < 


(log ir 


and 


L 


■a„+5 cs 

sin ipdx\ < 

-8 


log | 


for a lifting ip of m and for n = 1,..., TV, where C = C(a , TV, R, Co)- The first inequality is an 
immediate consequence of ©ID and ©3D- and the second follows with the same arguments as in the 
proof of Lemma 1131 These estimates are similar to (1421) and (1431) . but now we know that they hold 
true for non-minimizing configurations (under the energy control (1611 ) ) and the Pohozaev identity 
previously used for the pivof of ©2D and (1551) is no longer needed. 
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Remark 31. Inequalities (1611) and (1631) further imply that 


[ \S7u\ 2 dx< 2Co + ( ^ 1 . 

Jb S M ~ (log I ) 2 

If we combine this estimate with m, then we also obtain 

f |Vu| 2 dx < 1 (2C 0 + 3Ci + 2 [ |V< d \ 2 dx j 

JB*(a) (logf) V JB*{a)\Bt(a) ) 

for any r > 0. Furthermore, since u* a d is known quite explicitly, the last integral is typically not too 
difficult to estimate. 


6 Proof of the main result 

We now prove Theorem [2] To this end, fix a £ Am and d £ {±1}^. Set 7 n = d n — cos a for 
n = 1,..., N and 

N 

r = E^- 

n—1 

Furthermore, set 7± = l=Fcosa and recall Definition l26l which introduces the function e : {±1} —> R 
with 

e(±1) = l7o(( log s) 

We divide the identity from Theorem [2] into two inequalities, which are proved in Section [6. 2 1 and 
Section 16.31 respectively, after some preparation. Throughout the proof, we indiscriminately write 
C for various positive constants that depend only on a, N, a, d, and occasionally on the exponents 
of L p -spaces appearing in the context (always denoted by p or q). 


6.1 Preparation 

Define 


for x £ R+. Recall the functions 


wq(x) = arctan 


X2 \ 7rai 


xij 2 |xi| 


N 


N 


U *a,d = E &nd = E 

n=1 n =1 

from Section [2j Consider a number r £ (0, p(a)]. For n = 1,..., N, let 

An = In log(2 - 2a 2 ) + ^2 lkPa k (a n ) 
k^n 

again and recall estimate (l36l) . which implies that 


H*aA X l) - An - 7 nlog- 


\Xi CLr\ 


< Cr 


(71) 
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for x\ £ \a n — r, a n + r]. Also define 

Wn = ^7 ku ak (a n )- 

k^n 


Then 


and 


I u*a,d( x ) -w n - lnWo(xi - a n ,a,’ 2 )| < Cr in B r (a n , 0) 


(72) 


\\7u* ad (x) - 7„Vw 0 (^i - an,x 2 )\ < C in B r (a n ,0) (73) 

for 1 < n < N, because u* d {x) — oj n — 'y n wo{xi — a n ,X 2 ) is a smooth function that vanishes at 

{OjTI ) 0 ) • 

We now study how <5 changes when we replace e by e/r for a number r £ (e, 1], since we will have 
to rescale the magnetisation about the centres of the Neel walls. We have 


/ e r \ ( log r \ 

log I - log -J = log 6- logr + log II - I . 

Since log(l — £) < — £ for £ £ (0,1), we obtain 


log log - j < log 6 — log r — 


logr 

loge' 


Similarly, if e is sufficiently small (for a fixed r), then 


log log -^ > log 6 — log r 


2 logr 
loge ' 


(74) 


(75) 


6.2 A lower bound for the interaction energy 

The purpose of this section is to prove the inequality 

//f\2 -A n 

M Si nf ( (‘° g 6 ) E ‘ - ~ ‘° g 1 ) - g e(<W + W[a ’ d) (76) 

for the function W defined at the end of Section [5J which amounts to half of the statement of 
Theorem [5] 


First step: use minimisers Clearly it is sufficient to consider functions m e £ W 1,2 ((— 1,1); S 1 ) 
that minimise E e in M(a : d). Then 


by Proposition l27l 


lim sup 

e \,0 




(77) 


Second step: prove convergence away from the walls This part of the proof is similar to 
the proof of Theorem [22] Let tp e £ W 1,2 (— 1,1) such that m f = (cos</j e ,sinc^ e ). Define u f = U{m e ) 
and v e = u e log i. Then by Theorem [2H1 we have a sequence £fc \ 0 such that v ek — 1 v weakly in 
W’ 1,2 (fi r (a)) for every r > 0 for some function 

v£u* aid + W 1 ’ 2 (R 2 + ). 
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In fact, for any fixed r > 0, owing to Lemma [24] (with fl = R(j_), Theorem l28l and Remark l30l there 
exists a number /3 > 0 (depending on r) such that |sin</? e | > /3 at distance at least r away from 
— l,ai,... ,ajvj 1 for e small enough. For r, R > 0, define X ri ft(a) = (f l r (a) D B^{0))\(B r {—l,0) U 
-B r (l,0)). Then we can use Lemma |1 II and standard elliptic estimates to obtain uniform estimates 
in IF 2,2 (X rj fl(a)) for any r,R> 0 and any sufficiently small e. Therefore, we even have v ek v 
weakly in W 2 ’ 2 {T, r ^n(a)) for all r, R > 0. Furthermore, we have v tk —> v uniformly in R+\.E?2(0) by 
standard estimates for the Laplace equation. It follows that limi x i_ >co v(x) = 0. 

Obviously Av = 0 in and =0 on (—1, — oo) x {0} and on (1, oo) x {0}. By Lemma [TT1 
we also have 


lim sup 

e\0 


log 


0 •/. 


an+i-r 


(</?") 2 dxi < oo 


,+r 


for n = 1 ,,N — 1 and any r > 0, and we have similar inequalities in (r — 1, a\ — r) and in 
(ajv + r, 1 — r). Since 

u'e = e<p"/ sin(p e 

by (1371) . we conclude that u(-, 0) is locally constant in (—1, l)\{ai,..., a at}. But there is only one 
function in the space u* a d + VF 1,2 (R^_) with these properties (which can be seen with the arguments 
from the proof of Theorem 1771) . and thus we have 


u — u a,d 111 "“■+' 

Since |^(-, 0) = — m' u log \ on (—1, l)\{ai,..., a at}, it also follows that there exists a sequence 


dv, 
dx_ 

efc \ 0 such that 


(mu k - cos a) log -- > v 

Ok 


locally uniformly in (— 1 , l)\{ai, ..., a at} for a function v : (— 1 , 1 ) — > [— oo, oo] such that /i* d — v is 
locally constant in (—1, l)\{oi,..., oat}, where Sk = £fc log With the same arguments as in the 
proof of Theorem [27J we show that v = fi* . and 


(mi e — cos a) log — —> /u* d locally uniformly in (—1, l)\{ai,..., a at}- 
o ’ 


Now for any r £ (0, p(a)], we have 


I Vu 


a.d I 


' dx < lim inf 

/n r (a) 

Furthermore, by (ED, we have 


[ |Vw e 
Jn r (a) 


I 2 dx = lim inf 

e\,0 


>°4 


/ Q r (a) 


|Vw e | 2 dx 


mi e (a n ± r) — cos a — 


7n log y + A„ 


log; 


< 


Cr + o(l) 


log t log 


1 • 


(78) 


(79) 


Here and subsequently, we use the notation o(l) for any quantity that converges to 0 as e \ 0, with 
a rate of convergence possibly depending on r. 


Third step: rescale the cores Fix n £ (1, and r £ (0,p(a)], and define the functions 

m e : (—1,1) —> S 1 and u e : Ri —> R and the number e by 

m e (x i) = m e (rxi + a n ), 

Ue(x) = U e {rX 1 + a„, VXz), 
e 

e = -. 
r 
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Then we have 


and 


to'I dx i 


efjrn'fdx,=e£ 

I \X7u e \ 2 dx= I \S7u t \ 2 dx. 

Jb+( 0) JB+(a n ,0 ) 


(80) 


(81) 


Moreover, if we denote u*(x ) = u* d {rx\ + a n ,rx 2 ) and v £ = u e log then by the observations in 
the second step, we have v e — u* —*• 0 weakly in W 2 ’ 2 (B^ (0)\B 1 / 2 (0)). In particular, if we fix a 
number q > 2, then we have strong convergence of the boundary data in W 1,q {d + B+ (0)). Because 
of (1731) . we have 

\\x ■ Vr(a0||z,~ ( 3+B l( o)) < Cr. 

Hence 

Cr + o(l) 

IF ' VM e||L9(a+Bi(0)) < --• 

From Theorem [75] Remark [31] and inequality (1731) . we also obtain the inequality 


V« e — 


7rF 


M 2 log 7 

where 5 = e log i. We then also obtain 


|Vfi e || 


< 


c 


L 2 (B+(0)\B 5 -(0)) 


L*(B+m ^ log i» 


\x ■ Vu £ ||i<i(9+ Bl (0)) < 


Cr + o(l) 

logi 


and 


VfL — 




|a:| 2 logi 


I Vu £ || 


L*(B+m < , 1 • 


c 


L2(B+(0)\B- S (0)) 


log- 


Because of this and urn we may apply Corollary [21] to d n rhi e with 


7 = dn7n, r] = Cr + o( 1), and C = d n ( X„ + 7„ log -]+ Cr + o(l). 


In view of Definition [70] we conclude that 



> 7T7n l°g 7 + 2 e(d n ) - 27T7^ log - - 2n^ n X n - Cr - o(l). (82) 

S r 

We recall (1771) and (1751) . Now (1501) . (15T1) . and (1571) imply that 



> 2e(d n ) - 2ir^ 2 log - - 2n^ n X n - Cr - o(l). 
r 


46 
















Since this means in particular that 


l f r a n+r 

log i ' 


m'\ 2 dt + 


I B+(a n ,0) 

and since we have dzzD, the inequality also yields 


|Vm £ | 2 dx > Try 2 +o(l), 




s 

J a. 


a n +7* 


\m'\ 2 dt 


/s+(a n , o) 


|Vu e |“ dx - Try 2 log - - log - 


1 


1 


S r , 

> 2 e(d n ) - 2TTj n X n - Cr - o(l). 


Fourth step: combine the estimates If we sum over n and use (1751) . we obtain 

p \ 2 / p p \ N p n 

log ^j E e (m e ) - (log- -log- ) > J2(e(d n ) - 7r 7 „A„) + - / |V< jd | 2 dx - Cr - o(l). 

/ \ / n= 1 J 


Thus 


( / p \ 2 ttT 1 i ^ p r 7tT 1 

^log-J E e (m e ) — log — J > J^(e(d n )-7Tj n A n ) + - J^ ^ ^\S7u* ad \ 2 dx-^-\og--Cr. 


Finally we let r \ 0. Recall that 
1 


lim 

2 r\,0 


( [ \Wul d \ 2 dx-nT log- ) = E* d (u* d ) = Wi(a, d) 
\Jn r (a) r J 


and 


N 


-^'^2,lnX n = W 2 (a,d) 


for the functions W± and W 2 defined in Section 12.31 Hence we conclude that 

limmf ^ ^log ^ E e (m e ) - log ^ e(d n ) + W (a, d). 

That is, inequality (1761) is indeed satisfied. 

6.3 An upper bound for the interaction energy 

We now want to prove the inequality 

/ / f\ 2 p y \ N 

limsup (log - ) inf E e - log - < V e(d n ) + W(a, d), 

e\0 yv °) M{a,d) 2 dj ^ 

which complements GSD. 


(83) 
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First step: glue energy minimising cores into the tail profile Define 


r A n + 7 nlog i 

K n = -1-1-' 

In log £ 

Fix r G (0, p(a)\ small enough so that n r n G (0,1) for sufficiently small values of e. Choose minimisers 
G M| 7n | of the functionals ui 7 ”' and let w" be the solutions of the corresponding boundary value 
problem (l47l) (l49l) with 

[ u n t dx = 0. 

Vb+(0)\b 1/2 (0) 

Dehne /z" = d n \x ” and zt” = 

Let zy G C 00 (K) with 77 = 1 in (—00, |] and rj = 0 in [|, 00). Set 


<Pn(xi) = V 


\x\ - a n 


and ip n {x) = ?7 


\J{x 1 ^ an) 2 + : 


Note that = 0 on R x {0}. Now define 


mu{xi) = cosa + 


d*g,d( X l) 

log i 


+ Pnfcl) f(l - «£)/Xe/r ( 

n=l A V 


Xi - a r 


+ K^d n — cos a — 


K,d( x 1) 
iogy 


and 


0*0 = + J2 Mx) f(l- <K/ r ( 

® 5 n—1 ' ' 


(xi - a n ,x 2 ) 


a,d 


(x)~‘ 


!°g5 


Then m\ t (a n ) = d n for n = 1,..., N. If e is sufficiently small, we have —1 < mi e < 1. Hence there 
exists a function m 2e such that m e = (mi e ,m 2e ) G M(a,d). Let u e = U{m e ). 


Second step: estimate the magnetostatic energy in terms of u e Since u* d is harmonic on 
K+, we compute 


Au e (x) = A ^n(x) f( 1 - <)< /r (- 

n—1 ' ' 


(xi -a n ,x 2 )\ K,d( x ) ~ Ur 


log I 


,oY^V 7 / ( \ ( —1/1 r m n fi x l-a n ,X 2 )\ ^ U l,d( X )\ 

+ 2^Vz Pn(x)^r (1 - K n )\7u e/r ^--- j ~ lQg 1 j ■ 

Let = B^ r ^ 4 (a n ,0)\B r / 2 (a n ,0). Using Theorem [22] and the inequalities (1721) and (|73|) . we infer 


l e/r 


(x 1 -a n ,x 2 )\ K,d( x ) ~ 


log 


< 


and 


Vzt' 


e/r 


(x\ — a n ,x 2 )\ rVu* a4 (x) 


log I 


Cr + o(l) 
L °°(^n) ^g 5 

Cr 1+2 /P + o(l) 


< 


LP(E^) 


log 7 


(84) 


(85) 
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for any p < oo and all n = 1,... ,N. (Here o(l) again stands for any quantity that converges to 0 
as e \ 0, with a rate of convergence that may possibly depend on r.) Therefore, we have 

„ . _ „ Cr 2 ^- 1 + o(l) 

l|Au e LP(R2) < - : -J-• 

+ logy 

We also compute 


du e dx 1 ^a,d( X ^) A 

— (071,0) =--- 1 - >Vn{ x l) 

dx ^ lo g? (S 

N 


d ^/r ( 

( xi - a n \ 


\ r dx 1 

K r ) 

log| J 


= -m' u {xi) + ^2p' n (x i) ("(1 - <)/i"/ r ( 

n=l ' ' 


xi - a r 


+ dniil, — cos a — 


dajxi) 
1°S 7 


Using Theorem \22\ we then see that 


//"(xi) — cos a + 


7„log|xi| 


lo g? 


< 


o(l) 


!°ga 


for x\ G [— §, — g] U [3, f]. This, together with (f7T]h implies that 


(! - O+r 


Xi - a r 


+ d r K r „ — cos a — 


V*a,d( X l) 


lo g 7 


< 


Cr + o(l) 

log* 


for any aq £ [a„ - a n - §] U [a n + §, a n + ^]. 
Recall that u f is the solution of 


*■+’ 


A u e = 0 

1 

d£ = -m'ie on (—1,1) x { 0 }, 


du e 

&X2 


Thus we have 


= 0 on (—00,1) x {0} and (1, 00) x {0}. 

C7.2/P-1 + 0 (i) 


||A(u e -fo)|| iP(R p < 


for an arbitrary (but fixed) p £ (1,2) and 

d 


8X2 


(' U t ~ Ue) 


< 


log | 

(7 + 0(1) 

log? 


( 86 ) 


(87) 


( 88 ) 


Also note that the support of A (u e — u e ) is contained in B*(a ) and the support of -(u e — u e ) is 

contained in {J^ =1 {a n — r, a n +r). Thus if A4(R 2 ) denotes the space of Radon measures on R 2 , then 
after extending to R 2 by an even reflection on R x {0}, we have 


whence 


II A(u e - Ue) || A+R 2 ) 


Cr + o(l) 
lo g? 


||V(u £ Ue)|| i9 ( B +( 0 )) A 


Cr + o(l) 

lo g ? 
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for any fixed q £ [1,2). Theorem QT] and (l84l) imply that 

C + o{ 1) 


I || / al n+rn" < ^ 




It then follows that 


|VU £ || 


L«(B+( 0)) ^ 


logl 

C + o(l) 
lo §7 


(89) 


: dX. 


as well. We will use this inequality for q = in conjunction with (1871) . 
Let 

Then it follows that 

f |Vu e | 2 dx = — f 
Jr 2 , J( 


l 


B+(0) 


/R2 


/ 

Jr * 


(—1,1) x {0} 

Vi • V« f dx — 


, _ x du e 

( U e - U e )~ - dx 1 

0X2 


f (u e — u e )— — (u e — u e ) dx 1+ f (u e — u e )Au e dx. 

J ( — 1,1) X { 0 } OX 2 J R 2 


Now with the help of 
derive the estimate 


, and the continuous embedding 0)) L q ^ 2 q \— 1,1), we 

f t - ^ d ( - ^ j ^ C(r 2 ~ 2 / q + o(l)) 

/ (U e -U e )- (U e -U e )dx 1 < - ; - —5 -■ 

'(- 1,1)X{0} ^2 (fog l) 2 

(Note that for q = :i ' 2p _ 2 , we have 2 — | = | — 1.) Furthermore, by (l87l) . (|89l) . and the Sobolev 

C(r 2 / p - x + o( 1 )) 


inequality, 


/ ( 


— u e )Au e dx < 


MY 


If we choose p = |, then we obtain 


[ \\7u e \ 2 dx< [ |Vu t | 2 dx+ °^ + 0 ^ 
Jr 2 Jr 2 (log i 


(log 7)^ 


Hence 


E e {m e ) < ^ f \m' e \ 2 dt+^- f \Vil e \ 2 dx+ C ^ + 0 ^ 
1 J-i J Jr 2 (log 4 




(90) 


Third step: estimate ||Vu e || i 2( R ^) We clearly have 

f |Vfo| 2 dx=—^f | V< 

Jn r (a) (log ^) Jn r (a) 

In 13+ ( a n , 0) for n = 1,..., N, we have 
Wu e (x)=r~ 1 ( l-<)V< /r . 

+ {^n{x) - 1) ^r _1 (l - <)V 

+ Vi n (x) ("(1 - <)< /r 


* d\ 2 dx. 


(91) 


(xi - a n ,x 2 ) 


(xi -a n ,X2)\ V u *a,d{x ) 

log 4 


(xi - a„,a;2)\ K,d( x ) ~ 


lo §7 
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Thus using (l84l) , (l85l) and Theorem [22l and observing that 


|V<d| 2 


< 


c 


6/rlli 2 (B+(0)\B 1/2 (0)) - 


( l0 g lY 


by Theorem 1171 we obtain 


II V “ £ IIl 2 (B+(«„, 0 )) - ll VM e/i'llL 2 ( B +( 0 )) + ^ og l \2 ■ 

Combining (1DT1) and (1921) . we now find that 

N 

X7u e \ 2 dx <-—^—2 f \\/u* ad \ 2 dx + J2{1 - 2n r n ) f |V< /r | 2 dx + 
(logy) Jn r (a) n=1 Jb+( o) 

f iv«; 

Jb+( 0) 


( 92 ) 


Cr + o(l) 

N |) 2 ' 


Since 


. e/r | ! d*>S 


n 12 


c 


logj (log±)' 


by Theorem [TT1 it follows that 


|V£L 


N 

'dx< 1 2 [ |V < d | 2 dx + J2 f |Vn" /r | 2 dx 

(logy) ■'fir (a) n=l 0) 


N 


_ ^ 27T7 n (A w + 7 nlog^) + Cr + o(l) ^ 


n —1 


(log 


( lo gy)‘ 


log - 


Fourth step: estimate the exchange energy Note that we have |mi e — cos o:| < 

(—1, ai — |], in [a n + %,a n +± — |] for n = 1,..., N, and in [a/v — 1) by Theorem \7I\ Moreover, 

by m, we have 

dn - mie(xi) = (1 - k£) - Me/r 

, , T 

— cos Of —--j-1- d n n r 

lo S 7 


X\ — a r 


+ (1 - tfiniXi)) ( (1 - K r n )l% /r 

x\ — a. 


xi - a r 


= \d n - ii 


and 


e/r 


d n + mie(zi) = (1 - O ( d n + n 


i-m 
log V 


e/r 


x\ - a r 


+ (1 - <p n (x i)) - l)/z" /r 

xi — a x x 


H - ‘ld n l < i. 
x\ — a. 


— dnkL + cos a + 


V*a,d( x l) 

log | 


— I d n + /i) 


/r 


1-^ 

i°gy/ 
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in (a„ — r, a n + r). We also have 


4e0i) = 


1 - KL d ^/r (XX- 


dx i 


+ <p' n (x i)((i-<K/ r ( 

- (1 - <p„(a:i)) 


Xi — a„ 


+ d n Kl, — cos a 


logl 7 


( 1 - < d ^/r 

/ad - a n ^ 


l r dx 1 

{ r ) 

iogl J 


near a n . We have 

du* du* j 

J r A (xi) = -^r A (xi,o)- 

CLX i <7X2 

Hence defining = ( a n — r,a n — |) U (a n + a n + r), we have 


1 - d^ /r 

( Xl - «n ^ 


r dxi 

l r J 

log I 


< 


L 2( T r) 


cy? + o(i) 

iogi 


by Theorem [55] and ([75]) . Recalling ([55]) , we then compute 


N 


-1 


<i 2 = 


rl (sr^"/r) 2 _ dt+ o(l) 


S 2r J -1 1 - (M^) 2 (log A) 2 ' 

Combining (l90l) . (l9fil) . and (|94]l . we now find 

w 7T7n(A„ +7„logi) C' v / f + o(l) 


(94) 


^ 1 !\ 2 / l V <dl 2 dx + Y, jnf d/r - IZ ■ 

2 (log t) •'tV(a) n=l lTnl n=l 




( lo g iY 


Fifth step: estimate the core energy Recalling Definition [26] we see that 

/r' /2 

e/r “ 2bg ■ 


inf E h r \ < ^Hn + e(dn) + o(l) 

A/f e/r — o 1 

M hn I 


( lo U) 


2 ’ 


where 


6 r 

«= -log-. 

r e 


Using the estimates ((731) and ([75]) . we obtain 


Mi 


i„f E'r I < + e(<i " ) + ° (1) . 

I™1 Er 2 log 3 2 (log 3) (log 3) 


Sixth step: combine the estimates It follows that 


2 N 

^log ^ E e (m e ) < i ^ |V< )d | 2 dx + ^- log i log i + yj e{d n ) + W 2 (a, d) 


+ C-v/r-4- o(l). 
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That is, 


lim sup 

e\jO 



¥ lo 4) - If i vu m! 2 dx - ~v iog \+ e ( d ™)+ w2 ( a ’ d ) 

/ j£l r {a) n=l 

+ Cy/r. 


Letting r \, 0, we finally obtain 


(( i \ 2 . ^ j \ v 

ihnsup f (^log - J E e (m e ) -iog - J <^e(d„) + I V(a,d), 

which amounts to inequality (l83l) . This completes the proof of Theorem [2] 
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